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Abstract 

High-energy factorization in QCD is investigated beyond leading order and its 
relationship to the factorization theorem of mass singularities is established to 
any collinear accuracy. Flavour non-singlet observables are shown to be regu- 
lar at small x order by order in perturbation theory. In the singlet sector, we 
derive the relevant master equations for the space-like evolution of gluons and 
quarks. Their solution enables us to sum next-to-leading corrections to the 
small- a; behaviour of quark anomalous dimensions and deep inelastic scatter- 
ing coefficient functions. We present results in both MS and DIS factorization 
schemes. 
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1. Introduction 

Hadronic processes at large transferred momentum p t are accurately investigated by 
using perturbative QCD. The comparison between theoretical predictions and jet physics 
data from high-energy colliders has enhanced our confidence in perturbative QCD up 
to the 10% accuracy level [|J. The reason for this success is that not only the non- 
perturbative (higher- twist) contributions vanish as powers of A./p t in the hard-scattering 
regime pt ^> A (A being the QCD scale), but mainly the fact that logarithmic correc- 
tions to the naive parton model (i.e. lowest-order perturbative QCD) are systematically 
computable and, in most cases, known as a power series expansion in the 'small' (due to 
asymptotic freedom) running coupling a,s{Pt) ~ (A) In Pt/A 2 ) -1 . 

Precise quantitative tests of QCD and searches for new physics at present and future 
hadron colliders are nonetheless carried out at an increasingly large centre-of-mass energy 
\fS, thus opening up a new kinematic region characterized by small values of the ratio 
x = pf/S (S ^> pi). In this small-x regime, our capability to make perturbative QCD 
predictions decreases. The reason is twofold. First, the parton densities f(x,pf) of the in- 
coming hadrons (which are necessary as inputs in any perturbative calculation) are poorly 
known at very small values of x. Second, the perturbative expansion is slowly (or badly) 
convergent because of the presence of logarithmic corrections of the type ag(pf) ln m x: 
higher powers of atsiPt) associated with multiple hard-jet production can indeed be com- 
pensated by large enhancing factors ha.(S/p^). Note that this second point also affects the 
first one. In the case that a reliable set of small- a; parton densities can be extracted from 
a single experiment, we should be able to predict accurately their perturbative evolution 
with p\. 

The discussion above makes clear that our understanding of QCD in the small- X regime 
has still to be improved. In particular, on the theoretical side, it is highly desirable to 
keep under control and to evaluate reliably the QCD component which can be computed 
perturbatively. This can be done by combining the customary perturbative approach to 
hard-scattering processes with an improved perturbative expansion which systematically 
sums classes (leading, next-to-leading, and so on) of small- x logarithmic contributions to 
all orders in a $■ 

At present, the QCD multiparton matrix elements have been computed to double 
logarithmic accuracy in the small-a; region H||. They can be used to study the structure 
of hadronic final states in small- x processes, thus predicting some new distinctive features 
[0-^] such as the increase of particle (jet) multiplicity and the suppression of large rapidity 
gaps. Only some phenomenological investigations have been carried out 0, and more 
detailed analyses are needed. In this paper we do not consider the issue of the structure 
of the final states but rather concentrate on the evaluation of higher-order corrections for 
total cross sections. 

The leading high-energy contributions to total cross sections are single-logarithmic 
terms (aslnx)™ (higher powers of In a; cancel in this case) due to multiple gluon exchanges 
in the t-channel. In the case of hard processes which are directly coupled to gluons in 
the naive parton model, these leading logarithmic contributions can be resummed to all 
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orders in perturbation theory by using the high-energy or fc^-factorization theorem JS-TU]. 
The basic idea is to replace the collinear (or parton pole) factorization by gluon Regge 
pole factorization. 

Considering, for instance, the simplest small- a; process initiated by hard-gluon inter- 
actions at the Born level, namely the heavy-flavour photoproduction process 

7(p0 + Hp) - Q{pb) + Q(pd + x , (l.i) 

one can write the total cross section in the following factorized form || 

4M 2 a(x, M 2 ) = [ d 2 k f 1 — a(x/z, k 2 /M 2 , a s (M 2 )) T{z, k) . (1.2) 

J Jx Z 

Here the heavy flavour mass M (M ^> A) defines the hard scale of the process and 
x = AM 2 /S (S ~ 2 Pl • p > M 2 ). 



In Eq. (|1.2f) a is the basic high-energy hard cross section for the subprocess / y + g(k) — > 
QQ, computed to the lowest order in as a function of the transverse momentum k of 
the incoming off-shell (essentially transverse k ~ zp + fc 2 ~ — k 2 ) gluon g(k). On the 
other hand, !F(z, k) is the unintegrated gluon density of the incoming hadron h(p) and is 
related to the customary gluon density f g (x,M 2 ) via /^-integration 



r M- 

xf g (x,M 2 )~ / d 2 k F(x,k) . (1.3) 

J 

Therefore the /^-dependent factorization in Eq. ( |1.2| ) reduces to the leading-order collinear 
factorization |TTJ for S ^> M 2 ^> k\. However it holds also for S ^> k\ ~ M 2 , thus 
controlling all the logarithmically-enhanced terms (aslnx)™ associated with hard-gluon 
radiation in the final state. 

The resummation of the leading In ^-contributions follows from noticing that the hard 
cross section a is well-behaved at high energy (i.e. a(x, k /M 2 ,as) ~ x modulo lnx- 
terms, for x — > 0). Therefore the large perturbative corrections (aslnx)™ in the cross 
section ( |1.2|) are generated precisely by the /^-integration from the ones in the gluon 
density JF(x, k), as given by the Balitskii-Fadin-Kuraev-Lipatov (BFKL) equation ]T2" 

F(x,k) ~ -L e - Xlnx (k 2 )-^ , (1.4) 
A = 4C A — In 2 . (1.5) 

7T 

Inserting Eq. (|1.4| ) into Eq. (|1.2j ), one obtains the following perturbative result for the 
total cross section at very high energy || 

AM 2 a(x,M 2 ) ~ x~ x (M 2 )5 fc(l/2) , (1.6) 
, . x 1 /■«> dk 2 ( k 2 \ * /-i , 9 , , , , . 
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The main features of Eqs. ( |1 . 7| ) , derived from the fc^-factorization formula ( |1.2| ), 

are the following. The total cross section increases at high energy with a universal 
(process-independent) power behaviour S x , A being the intercept of the perturbative 
QCD pomeron in Eq. ( |1.5| ). This result is the consequence of the very steep behaviour 
(|1.4j) of the gluon density^] at small x and large k±. The normalization of the total cross 
section instead depends on the process and the process-dependent factor h(l/2) derives 
from the detailed and calculable transverse-momentum dynamics of the hard subprocess. 

The resummation of the leading-order contributions at high energy is a crude (although 
mandatory) approximation. This is true from both the theoretical (the perturbative QCD 
pomeron violates unitarity to leading order) and phenomenological (terms of relative 
order as are systematically neglected) sides. The evaluation of subleading contributions 
is therefore relevant i) to include corrections necessary to restore unitarity at asymptotic 
energies and ii) to estimate the accuracy and set the limits of applicability of the leading- 
order formalism. 



Unitarization effects have been extensively studied in Refs. [pr3Hl6[|. Although a system- 
atic calculational approach based on first principles is still missing, the likely conclusion 
emerging from these studies is that the full restoration of unitarity can be achieved only 
after the inclusion of higher-twist corrections. 

As regards the leading-twist contributions to hard-scattering processes, a calculational 
program of high-energy logarithms based on Regge behaviour is being pursued by Fadin 
and Lipatov ||17|| , and the evaluation of the two- loop correction to the BFKL kernel now 
seems feasible. 

In this paper we follow a different approach towards the computation of subleading 
corrections at high energy. We show how the fc^-factorization theorem can be extended 
beyond leading order in a consistent way with all-order (leading-twist) collinear factoriza- 
tion. This allows us to set up a systematic logarithmic expansion both for hard coefficient 
functions and parton anomalous dimensions. Moreover, once these quantities have been 
computed to a certain logarithmic accuracy, they can unambiguously be supplemented 
with non-logarithmic (finite-x) contributions exactly calculable to any fixed order in per- 
turbation theory. Note, also, that a further advantage of this approach is that the effects 
of the running coupling can be included exactly (at least in principle), thus avoiding the 
infrared instabilities encountered in phenomenological attempts P~5|JT~5|JT9[ ] to extend the 



BFKL equation beyond leading order. Obviously, with the present attitude, we abandon 
any demand to predict the absolute behaviour of the cross sections at very high energies 
(very small x) because higher-twist corrections are systematically neglected. However, 
since logarithmic scaling violations are systematically under control, we think that such 
an approach can be useful in making quantitative phenomenological predictions at high 
(but finite) energies and large transferred momenta 



The outline of the paper is as follows. In Sect. 2 we first recall the general framework 
of leading-twist collinear factorization. Then, on the basis of power counting arguments, 



*Note, however, that if the gluon density has a non-perturbative component steeper than the pertur- 
bative one, the former dominates over the latter at high energy. 



-4- 



we show how the high-energy factorization of Ref. || can be extended beyond leading or- 
der and consistently matched with all-order collinear factorization, in terms of resummed 
anomalous dimensions and coefficient functions. In Sect. 3 we start our calculational pro- 
gram in dimensional regularization by deriving the master equation for the high-energy 
behaviour of the gluon forward scattering amplitudes and computing the ensuing anoma- 
lous dimensions and normalization factor in the (modified) minimal subtraction scheme. 
The analogous calculation for the quark channel is performed in Sect. 4. Here, we obtain 
an algebraic equation for the quark anomalous dimensions to next-to-leading logarith- 
mic accuracy as(a,g In x) n , and present its explicit solution up to the six-loop order. In 
Sect. 5 we turn our attention to the calculation of the coefficient functions for deep inelas- 
tic lepton-hadron scattering. We compute both the longitudinal and transverse coefficient 
functions by resumming the logarithmic contributions ocsi&s In %) n - We also consider the 
all-order generalization of the DIS factorization scheme and obtain explicit resummed 
expressions for the corresponding quark anomalous dimensions with next-to-leading log- 
arithmic accuracy. Section 6 is devoted to summarizing our main results, and may also 
serve as a guide for the reader mostly interested in extracting perturbative QCD results for 
phenomenological applications. The definition of singlet and non-singlet parton densities 
and a few mathematical details are left to appendices A to C. 



Some of the results derived in this paper have already been presented in Refs. pT,22|. 



2. QCD factorizations 

This section is devoted to setting up the formal basis of the high-energy factorization. 
We start in Subsect. 2.1 by introducing the factorization theorem of collinear singular- 
ities and defining the parton densities and coefficient functions. Then in Subsect. 2.2 
we recall the proof of this factorization theorem to leading-twist order in the context of 
dimensional regularization. This formal apparatus is used in Subsect. 2.3 to develop a 
simple power counting at high energy and to show how high-energy factorization can be 
carried out consistently with all-order collinear factorization. The high-energy factoriza- 
tion formulae are fcj_-dependent and can be considered as the generalization of collinear 
factorization in terms of unintegrated parton densities and off-shell coefficient functions. 



Alternatively, high-energy factorization can be compared with Regge factorization p3 
The gluon Green functions and the two-gluon irreducible kernels entering as building 
blocks in the high-energy factorization discussed in Subsect. 2.3 are the analogue of the 
reggeon trajectory and its residue. However, this analogy has to be taken with caution 
since it does not properly account for the issue of factorization of collinear singularities. 
Finally, in Subsect. 2.4, we show how the high-energy factorization leads to resummed 
anomalous dimensions and coefficient functions. 



2.1 Hard processes and parton densities 



The factorization theorem of collinear (mass) singularities p4] , ^5| , p!T|| states that, in a 
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general hard collision (i.e. a scattering process involving a large transferred momentum 
< p 2 >= Q 2 ^> A 2 ) of incoming hadrons, all long-distance (non-perturbative) effects 
can be factorized into universal (process-independent) parton densities thus leading to a 
perturbatively calculable dependence on the hard scattering scale Q 2 . 

Considering, for the sake of simplicity, the case of a single incoming hadron (like the 
heavy-flavour photoproduction process (|1.1| ) or deep inelastic lepton-hadron scattering), 
one can write the dimensionless cross section F(x, Q 2 ) ~ Q 2 a(x, Q 2 ) as follows (x = Q 2 / S) 

F(x,Q 2 ) = J2 I' dzC a (x/z;a s ^ 2 F ),Q 2 /fi 2 F ) f a (z,fi 2 F ) . (2.1) 

a Jx 

Here C a are the pro cess- dependent coefficient functions, f a are the parton densities (a = 
qi,(ji,g, i = 1, • • • , Nf, Nj being the number of flavours) of the incoming hadron and fi F 
is an arbitrary factorization scale such that \j? F ^> A 2 . The observable F is independent of 
fi 2 F . Correspondingly the /independence of C a on the r.h.s. of Eq. (|2.1| ) is exactly cancelled 
by that of f a . Moreover, both the a^-dependence of the coefficient functions and the scale 
dependence of the parton densities are computable as a power series expansion in In 
particular, the parton densities fulfil the renormalization group evolution equations 

df a (x,/j 2 ) x ^ f 1 dz 



dln/i 2 



£/ -P ab (a s (v 2 ),z)f b (x/z,v 2 ), (2.2) 



where P a b(as, z) are generalized Altarelli-Parisi splitting functions which are computable 
in QCD perturbation theory: 



P ab (a s ,x) = j:[^) Pr'\x). (2.3) 



n=l 

The corresponding expansion for the coefficient function is 



C a (x ]as ,Q 2 /n 2 ) = (a s y 



n=i V2vr/ 



(2.4) 



where the integer power p depends on the process. 



The x-dependence of the factorization formula (|2.1| ) and of the evolution equations 
]2|) can be diagonalized by introducing the iV-moments of the cross section 

F N {Q 2 )= [ dx x N ~ 1 F(x, Q 2 ) (2.5) 
Jo 

and the analogous moments of any other function of x. It is also convenient to define the 
rescaled parton densities f a 

fa(x,{l 2 ) = Xf a (x,/J, 2 ) , (2.6) 

and the anomalous dimension matrix [] 



1ab,N 



(a s ) = [ dx x N P ab (a s ,x) = P ab . N+1 (a s ) . (2.7) 
Jo 



"Out definition differs from the standard one in which 7jv = Pj 



N- 
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In the iV-moment space, Eqs. ( |2.1| ),( |2~2] ) respectively read as follows 

F N (Q 2 ) = E Ca,7v(«5(/i|),Q 2 //iF) iU(/4) , 



din [J? 



J2lab,N{0is{^ 2 )) fb,N{^ 



(2.8) 
(2.9) 



Eq. ( [2.1|) is correct apart from terms vanishing as inverse powers of Q 2 at high Q 2 
(higher-twist corrections). In this regime, a consistent perturbative use of the collinear fac- 
torization formula Q2.1| ) requires the knowledge of {P^\ C^} (leading order), {P$, C^} 
(next-to-leading order) and so on. Therefore, QCD predictions for the cross section are 
usually obtained by computing in fixed-order perturbation theory both the splitting func- 
tions and the coefficient functions. The splitting functions up to two-loop accuracy have 
been known for a long time [^5| - |28|] . Next-to- leading order coefficient functions have been 
computed for most processes [29], and, in the cases of deep inelastic lepton-hadron scat- 
tering (DIS) |30| and Drell-Yan process the next-to-next-to-leading terms are 
also known. 

As discussed in Sect. 1, higher-order contributions to the cross section and, hence, to 
splitting and coefficient functions are logarithmically enhanced at small x. More precisely, 
in the small- x limit we have [] 



3>(n-l), 
ab 1 



1 r 



X) ~ 



X 



In™" 1 x + £>(ln 



n-2 



X 



X) ~ 



In™- 1 x + C(ln 



n-2 



X) 



lab,N{<ys) = E 
k=l 



—) A ab +ocs^—) B ab +... 



(2.10) 
(2.11) 



This small- x behaviour corresponds, in iV-moment space, to singularities for N — > in 
the form 



X12) 



C a ,N(a s ,Q 2 7 'fJ?) 



n P Mo) 

a 8 u a,iV=0 



) £?>(tfV) + --- 



1 + E(f) D a\Q 2 /V 2 ) 



(2.13) 



These singularities may spoil the convergence of the perturbative expansions ( |2.3|) and 
( |2.4| ) at small x (high energy). Nonetheless one can consider an improved perturbative 
expansion obtained by resumming the leading (A { ^,D^), next-to-leading (B§\eW), 
etc., coefficients in the high-energy regime. Once these coefficients are known, they can 
be combined with Eqs. ( |2.3D and ( |2.4| ) (after subtracting the resummed logarithmic 
terms in order to avoid double counting), to obtain a prediction throughout the region of x 
where as ln(l/x)^l (or ots/N&l), which is much larger than the domain asln(l/:r) <C 1 
where the a^-perturbative expansions fl2.3|) , ( |2.4j ) are applicable. 

^ According to our definition, the lowest-order coefficient function is normalized in such a way 
that C(°'(x) ~ x, modulo In a; corrections, at small x. 
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As outlined in Sect. 1 and discussed in detail in Ref. || (see also Sect. 5), the k±- 
factorization formula ( |1.2|) allows one to resum the leading high-energy contributions 
(aslnx)™ ((as/N) n ) both for the anomalous dimensions and for the coefficient function. 
By comparing Eqs. ( |1.2| ) and ( [2.1| ) and considering the relation ( \1.3\) between unintegrated 
and full parton densities, we see that the fcj_-factorization is, in a sense, more general 



than the collinear factorization in Eq. ( |2.1| ). The latter is recovered after the transverse 
momentum integration in the former. This relationship can be made explicit in a simple 
way to leading order || . However, beyond leading order, a careful analysis of the \ow-k± 
integration region in (|1.2| ) is needed in order to disentangle non-perturbative and higher- 
twist effects from the perturbative ones. 

This issue can be investigated from a formal viewpoint by considering cross sections 
at parton level instead of hadronic cross sections. In the partonic cross section, the 
non-perturbative contributions show up in terms of collinear divergences. Once these 
divergences are properly regularized, they can be factorized and subtracted by a procedure 
of renormalization of the bare parton densities. This procedure can be carried out to any 
collinear accuracy, although the ensuing anomalous dimensions and coefficient functions 
are no longer separately regularization/factorization scheme independent beyond one-loop 
level. 

The point we want to address in the following is that the fcj_-factorization can be 
consistently implemented beyond leading order without spoiling the all-order factorization 
of collinear singularities. In particular, the resummation of next-to-leading contributions 
at high energy can be performed by properly taking care of the factorization scheme 
dependence of the splitting and coefficient functions. To this end, let us first recall the 
formal basis of collinear factorization. 



2.2 Collinear factorization 

In order to present a formal derivation of the factorization theorem of collinear singu- 



larities, we follow the technique developed by Curci, Furmanski and Petronzio p5| . The 
dimensionless cross section F(x, Q 2 ) in eq. ( [2.1| ) is first expressed in terms of partonic 
cross sections and parton distributions in the form (Fig. 1) 

F = F<°>(. . . ,p) /W(p, . . .) + [F 4 (0) (. . . ; Pl ,p 2 ) fi°\ Pl , P2 ; ...) + ...] , (2.14) 

where, in the limiting case of on-shell partons (p 2 = 0), the first term on the r.h.s. picks 
out the leading-twist contribution we are interested in. Nonetheless, the partonic cross 
section is collinear divergent in the on-shell limit, so that both F^ and have to 
be regarded as properly regularized 'bare' quantities. 

We use the standard procedure of dimensional regularization, in which the bare cross 
section F^ is evaluated in n = 4 + 2 e space-time dimensions, considering (n — 2) he- 

^Here and in the following we use a symbolic notation in which the product C — AB of two kernels 
A and B understands the integration over the intermediate momenta (or the corresponding product in 
JV-moment space) and the sum over the intermediate parton species and their spin and colour indices. 
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licity states for gluons and 2 helicity states for quarks. The corresponding dimensional- 
regularization scale is denoted by /i. Once has been renormalized (using, for instance, 
the MS renormalization scheme) , collinear singularities are automatically regularized and 
show up as single poles in 1/e. The factorization theorem allows one to subtract these 
poles from and factorize them (to all orders in as) into process-independent transition 
functions T, according to 

F (0) = CT , (2.15) 

where the coefficient function C is finite for e — > 0. Using the transition functions V to 
define the 'physical' parton densities / 

/ = r/w, (2.16) 

one then recovers the factorization formula (|2.8|) by performing the limit e — > 0. 

The factorization procedure leading to eq. (|2.15|) is simplified if we evaluate the gauge- 
invariant partonic cross section F^ ' in a physical gauge. Denoting by p^ = P(l, 0, 1) the 
incoming parton momentum (Fig. 2), we introduce the following Sudakov parametrization 
for any other momentum k 

k » = zp » + k » ±+ k2 + k2 ? , ^ = (0,k,0) , p" = P(l,0,-l) , (2.17) 

and we choose the axial gauge p ■ A = 0, where the sum over the gluon helicities is given 
by the polarization tensor 

d^(k) = -<T -{• ——— P . (2.18) 
p ■ k 

Then we consider the expansion of F^ in terms of kernels and which are two- 
particle irreducible (2PI) in the t-channel: 

F (0) = C (0) (1 + K (0) + K (0) K (0) + ...) = C (0) g(o) ? (2 .19) 
g(0) = l + K (o) + K (0) K (0) + = _J__ {22Q) 
In the axial gauge in which we are working, the 2PI amplitudes are free from mass 



singularities [Q. Therefore all the collinear divergences originate from the integrations 
over the momenta coming out from the kernels and connecting them to each other in 
the (process-independent) bare Green function C/(°). The factorization formula ( j2.15|) can 
now be realized by introducing a suitable projection operator Vq which decouples G^> 
and in the spin indices and extracts the singular part of the d n k integrals (i.e. poles 
in e) thus decoupling and also in momentum space. 

For each kernel, one can write the decomposition K^ ' = (1 — Vc)K^ +VcK^°, where 
all the singularities are due to the second term on the r.h.s. Applying this procedure in 
an iterative way, one obtains 

g {0) = gr (2.21) 
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where all the e-poles have been subtracted from the 'renormalized' Green function Q: 
1 



1 - (1 - V C )K(°) 

and associated with the transition function Y 

1 



1 + (1 - V C )K® + (1 - V c )[K i0 \l - Vc)K^} + ... , (2.22) 



1 + (V C K) + (V C K)(V C K) + . . . 



1 - (VcK) 

k = k<® g . 

The coefficient function C in eq. fl2.15| ) is thus identified with 

c = c {0) g . 



(2.23) 
(2.24) 

(2.25) 



In order to show the consistency of collinear and fcj_-factorization (see Sect. 2.3), we 
have to recall the form of the projection operator Vc ||25||. Let us denote its action on 



helicity and momentum space respectively by and Pq, so that Vc = P( 



00 
C 



P^\ If 



A and B are two kernels connected by a parton of momentum k (Fig. 3), the action on 
the helicity space is 



AP^ B = A{...,k). 



.a'0' 



p 



.2 ) a >p \2p-kJ al3 
when the connecting parton is a quark and 



B a/3 ...(k, . . .) 



(2.26) 



A P { c s) B = A(..., — jtt (-gnB^.Xk, 



(2.27) 



when the connecting parton is a gluon. The operator P^) sets k^ = zp^ on the left- 
hand side (A), performs the dk 2 d n ~ 2 h integration up to the factorization scale \j? F on the 
right-hand side (B), and extracts the ensuing poles in e. 



Note from Eqs. (|2.26 ),( 2.27 ) that P^ acts on the left by performing the average over 



the parton helicities in n- dimensions. Note also that Pq is not (at least in principle) 
unambiguously defined. The factorization of Eqs. ( |2.15| ),( p.21| ) in terms of collinear- finite 
and divergent (for e — > 0) contributions can still be achieved if P^ extracts not only the 
£-poles but also any finite contribution for e — > 0. This leads to the factorization-scheme 
dependence of both anomalous dimensions and coefficient functions. The factorization 
scheme is completely specified once P^ has been uniquely defined. Equivalently, one can 
specify the explicit (and finite) e-dependence of the transition functions Y. One of the 
most commonly used scheme is the modified minimal subtraction (MS) scheme, in which 
the transition functions have the form 



Y ab , N (a>s(n F /n 2 ) e ,e) = P Q exp - / 

\e Jo 



ra s (n%/^rs s da 



— In [a 
a 



'2.28) 



In this equation we have reintroduced the explicit dependence on the parton labels and 
iV-moment indices. The symbol P a denotes, as usual, the path-ordered integration of 
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the anomalous dimension matrix 7 a fc,7v- Note the presence of the e- finite factor S e = 
exp{— £[^(1) +ln47r]} (ip(z) is the Euler ^-function), which characterizes the MS-scheme. 
The only formal simplification we have used in Eq. Q2.28|) is that we have considered the 
case of a fixed coupling constant as- As shown in the following, this simplification is 
sufficient for the purposes of the present paper. 



2. 3 Power counting and factorization at high energy 

The expansion in 2PI kernels introduced in the previous subsection is particularly 
useful to discuss the high-energy behaviour. High-energy (or small-x) logarithmic contri- 
butions are indeed generated by multiple gluon exchanges in the t-channel. Therefore we 
are led to consider kernels which are two-gluon irreducible (2GI). 

For the parton cross section (a = qi,q~i,g) , we single out the part which is 

2GI by selecting the first (starting from above in Fig. 4a) two-gluon intermediate state. 
Considering the small- N limit in iV-moment space, the 2GI kernel behaves as a p s (1 + 
as + ■ ■ ■), where the first term corresponds to the tree-approximation and the remaining 
terms stand for corrections which are subleading at high energy. The large perturbative 
contributions (as/N) k are thus generated precisely by fc-integration from the ones in the 
gluon Green functions Qjj£ (a = qi,q~i,g). In particular, since flavour non-singlet parton 
cross sections (App. A) get no contribution from pure-gluon intermediate states, we can 
immediately conclude that non-singlet anomalous dimensions and coefficient functions are 
regular for iV — > order by order in a$- 

A decomposition similar to that for can be performed also for the (flavour singlet) 
quark Green function (Fig. 4b). Since the 2GI kernel behaves in this case as as (l + as + 
. . .), we see that the quark anomalous dimensions contribute to next-to-leading terms 
a s( a s/N) k in the high-energy limit. 

Note that the expansion in 2GI kernels is more general than that in 2PI kernels. 
The 2GI kernels in Fig. 4a and Fig. 4b can still be expanded respectively as (C^ + 
Eb^gC^K^ + ...) and (K$ + Ect g Kg>KV> + ...). However, only the tree-level ap- 

proximations for and (b ^ g) contribute to leading order in as/N. A similar 
simplification does not occur for the gluon Green function Qj® because the gluon kernels 
Kf® contain terms of the type (as/N) k to any order k (modulo dynamical cancellations) 
in as- 

The expansion in 2GI kernels described so far allows one a simple power counting 
at high energy. The next step towards high-energy factorization consists of decoupling 
the 2GI kernels and the gluon Green functions with as/N fixed. This factorization is 
conceptually different from the collinear one, where, roughly speaking, one expands in as 
(or e) with as/s fixed. In particular, we cannot perform the collinear limit k 2 — > to any 
fixed order in as, because small- x contributions (as/N) h are associated with any value 
of k ||. On the other side, we do not want to spoil the collinear factorization, so that the 
high-energy factorization has to be valid for any value of e (i.e., in any number of space- 
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time dimensions). We are going to show that the fc^-factorization procedure introduced 
in Ref. || can be carried out consistently with all-order collinear factorization. 

The high-energy limit of the product A g Q(°J, A g being the 2GI kernel involved in the 
decomposition of Fig. 4, was discussed in detail in ||. Considering the decomposition of 
A g in Lorentz-invariant amplitudes, it was shown that the leading high-energy behavior 
can be extracted via fc^-factorization, i.e. by inserting into the two-gluon intermediate 
state a /cj_-dependent projection operator Vh as follows 

A 9 Gf}=A g V H Gf} + ... . (2.29) 

As in the case of the collinear projector Vc in Sect. 2.2, we introduce the notation Vh = 
Ph®Ph> where P^ and P^ denote respectively the action of the high-energy projector 

(s) 

Vh on helicity and momentum space. P H acts as follows 

A 9 P { h ] Gf} =A(..., k)f A -±^ {-g„)S^(k,p) , (2.30) 

whilst Pff sets k^ = zp^ + k^_ on the left-hand side (A g ) and integrates the right-hand 
side (Gj$) over the invariant mass k 2 at fixed k±. Note that the fcj_-dependence is left 
unaffected by Vh and, in particular, the 2GI kernel A g has to be evaluated with an 
incoming off-shell (essentially transverse k 2 = — k 2 ) gluon. 

Equation ( |2.29|) generalizes the fcj_-factorization formula ( |1.2| ) to n = 4 + 2 e space- 
time dimensions. The key point, however, is not just the formal resemblance between 
Eqs. ( |2.29| ),( pr3~0"D and (|TT2|), but rather the fact that Vh selects the correct high-energy 



behaviour in any number of dimensions. We mean that, for instance, in performing the 
approximation ( [2.29 ), we are not neglecting any contribution of order (as/N) h ■ e with 



respect to the four- dimensional case. This statement is a consequence of the fact that Vh 
is a 'true' projection operator: 

V H = V H , (2.31) 

and fulfils the property: 

Vh^Vc, (V h = V c iff k± = 0) . (2.32) 



Eqs. (|2.31| ) and ( p.32|) are self-evident for the momentum space components Pjj \ P@ and 
follows from the simple relations (<>0 denotes the average over the n — 3 azimuthal angles 
in the transverse momentum space) 

/ J^L i ^^ ± ki k± ^ dr{k = z P ) 

\-9»v)—£- = 1 , < —g- ><j> = ' ( 2 - 66 > 

for the spin components Pjj , Pc ' 

Eq. ( ^.32| ) guarantees the consistency between high-energy factorization and collinear 
factorization. Due to Eq. ( [2.32| ) we can first perform the high-energy approximation in 
Eq. ( p.29|) and then proceed to the all-order factorization of collinear singularities by 
applying iteratively the collinear projector Vc, as described in Sect. 2.2. 
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2.4 Resummation at high energy 

In order to describe how in practice this procedure works and leads to the high-energy 
resummation, let us consider the case involving a 2GI kernel which is collinear safe. In 
particular, we refer to the heavy-flavour photoproduction process ( |1.1| ) already introduced 
in Sect. 1. The corresponding 2GI kernel A g to lowest order in as is given in Fig. 5: the 
dashed line denotes the incoming on-shell photon and the full lines correspond to the 
heavy-flavour pair produced in the final state. Considering the case of an incoming-gluon 
partonic state and using the high-energy approximation in Eq. ( |2.29p , we immediately 
obtain the factorized formula (|1.2j ) (in n = 4 + 2 e dimensions) , which we can rewrite in 
the iV-moment space as follows 

AM 2 a N {M 2 ) = Jd 2+2s k d N {k 2 /M 2 ,as{M 2 /fy ] 6)^{k ] a s ,^,e)ff} ] ^,e). (2.34) 

Here a ~ A^(k) k^k^/k 2 , is the bare gluon distribution and, according to the 
action of the high-energy projector Vh, we have introduced the /c^-dependent gluon Green 
function 

r dk 2 

^(z,k;a s ,n,e) = J (-^£^(M) • (2-35) 



Equation ( 2.34| ) has to be regarded as the bare (and collinear regularized) version of 



the factorization formula ( |1.2| ), in the sense that it still contains collinear poles in e which, 
according to Eq. ( |2.21| ), can be factorized in the form 

T^\k; a s , fi, e) = ^gMots) ^ as] ^ £ ) r gg7N (a s {n F / n 2 ) £ , e) , (2.36) 

7TK 

where the MS-scheme gluon transition function is given by f\ 

( 1 rasSe da 1 
F gg ,N(as, e) = exp < - — 7 ra ,A r (a) } ■ (2.37) 



The function R N in Eq. ( |2.36| ) has no e-poles order by order in as, and cannot depend 
on the dimensional regularization scale \x in the limit e — > 0. Moreover, since the l.h.s. 
of Eq. ( p.36| ) is independent of the factorization scale /xp, the only k-dependence of Rjy 
allowed by dimensional arguments, for e = 0, is the following 

R N (k, fi F , a s ; n, e = 0) = R N (a s ) (k 2 /^)^^ . (2.38) 
The reduced cross section a in Eq. (|2.34|) is collinear safe because it corresponds to a 



2PI kernel of the type in Eq. (|2.19| ). Therefore, after using the transition function 
Tgg^ in ( [2.36 ) to 'renormalize' the bare gluon density as in Eq. ( [2.16 ), we can safely 



perform the e — > limit in Eq. ( |2.34f ) and obtain: 

4M 2 a N (M 2 ) = C N (a s , M 2 /fi 2 F ) f a , N (fi 2 F ) , (2.39) 

^We limit ourselves to the case of a fixed coupling constant as- In fact, the running coupling effects 
lead to subleading contributions as(as /N) k and can thus be neglected in the present leading-logarithmic 
analysis. 
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C N (a s , M 2 //i|) = h N ( lggiN (a s )) R N (a s ) (M 2 /^""*^ , (2.40) 

where the process-dependent part of the coefficient function Cn is given by the fol- 
lowing A;_|_-transform of the hard cross section a: 

dk 2 ( k 2 x 7 
/o "kMM 2 



h N { 1 )= 1 j —[ — ) a N (k 2 /M 2 ,a s ;e = 0) . (2.41) 



The result in Eq. ( [2.40|) gives the resummed expression (including the dependence on 
the factorization scale ii 2 F ) for the coefficient function to the leading order (as/N) k , 
provided the gluon anomalous dimensions 7 9g ,7v and the process- independent function R N 
are known to the same accuracy. A similar result was first derived in Ref. ||. The only 
difference with respect to [§] is that Eq. (|2.40|) takes into account the explicit depen- 
dence on the process-independent but factorization- scheme dependent function R^. This 
point is essential for precise phenomenological predictions at finite energies, when one 
has to combine the resummed coefficient function ( [2.40|) with fixed-order non-logarithmic 
contributions computed in a well defined factorization scheme of collinear singularities. 
Moreover this issue is relevant to extending the high-energy resummation to subleading 
orders, where a corresponding scheme dependence of the anomalous dimensions comes 
into play (see, for instance, Sect. 5). 

Note, however, that the presence of the process independent factor R N in Eq. ( |2.40| ) 
is no longer relevant (at least, to leading order) if one limits oneself to considering only 
ratios of cross sections. In this case R^ and the factorization scale dependence cancel 
in the ratio, thus leading to an absolute prediction in terms of hard scales and ^(7) 
functions of the type in Eq. ( |2.41| ) |Q . 



The high-energy contributions (cts/N) k are embodied in Eqs. (gjp , fl2T34T ) through the 
/^-integration of the gluon Green function Qfd- As discussed in and recalled in Sect. 1, 
in the case of four space-time dimensions the resummation of the leading terms (ces/N) k 
in the gluon density (anomalous dimensions) is accomplished by the BFKL equation []T2J . 
The analogous master equation in n = 4 + 2 e dimensions, which is necessary to compute 
both the anomalous dimensions and the function R^, is derived and discussed in the 
following section. 



3. The gluon Green functions 

The leading high-energy behaviour of the gluon Green functions Gg°J(k,p) and 
Gj®\k,p) can be easily derived by generalizing the soft-gluon insertion technique in 
Refs. PJ5|J34]] to the case of n dimensions. In the present paper we do not repeat all 
the detailed calculations described in f3|,|4j], but we simply sketch the main steps and 



properties which are necessary for the n-dimensional generalization. 



The starting observation is that the high-energy contributions (as/N) k to the gluon 
Green function are produced by radiation of soft gluons, that is real and virtual gluons 
carrying a very small fraction x of the longitudinal momentum p of the incoming parton. 
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Therefore, in the soft-gluon approach, the matrix element M^ k+1 \k,p), contributing to 
to the (k + 1)— loop order, is obtained from by the insertion of an additional 

(real or virtual) soft gluon with momentum q. Using the soft approximation for vertices 
and propagators, this insertion can in turn be factorized, leading to the recurrence relation: 

\M^(k,p)\ 2 = g 2 s {{M( k \k + q ,p)Vij£} t ( q )} 2 M( k \k + q,p) 

- [M«(A;,p)]t [j£7 t (g)] 2 M«(M} • (3-1) 

The explicit expressions for the real and virtual soft-gluon currents J^J t and J^l can 
be found in Refs. |||34[]. The point we want to emphasize here is that they do not explicitly 
depend on the number n of space-time dimensions in which the soft-gluon factorization 
is carried out. As a matter of fact, the n-dependence of the matrix element M^ k+1 ^ can 
only be due to the spin structure of the vertices (the scalar propagators i/(q 2 + ie) are 
the same in any number of dimensions). It turns out that a single (essentially eikonal) 
helicity flow dominates to leading-logarithmic accuracy in the high-energy limit. Thus, 
the relevant QCD matrix elements are the same as for n = 4 and the only difference 
comes from the n-dimensional phase space over which Eq. (|3.1| ) has to be integrated. 



Using the recurrence relation ( |3.1|) , performing the sum over the number k of loops 
and following exactly the same steps as in Refs. PJ3"1"|, one obtains an integral equation for 



the gluon Green function Q(°\ More precisely, considering the gluon density in Eq. ( |2.35| ) 
one finds (as = Cachs/^, Ca = N c being the number of colours) 



(27r/i) 2e 7rq 2 



• jVff^k - q; a s , //, e) - ^ J ^ ^(k; as, V, e)| • (3.2) 

Some comments are in order. In the case of n = 4 dimensions (e = 0), after azimuthal 
average over q, Eq. (|3.2|) reproduces the BFKL equation |12],[35|. Moreover, the integrand 



(not the phase space) of the homogeneous term in Eq. ( p.2[ ) is exactly the same as for 
the BFKL equation. As discussed above, this is a straightforward consequence of the dy- 
namical dominance of a single gluon polarization to the present accuracy. Note, however, 
an essential difference with respect to the BFKL equation: the full kernel of Eq. ( |3.2|) 
is not scale invariant. Indeed, scale invariance is broken by the dimensional regulariza- 
tion procedure and the eigenfunctions of the kernel are no longer simple powers (k 2 ) 7 . 
From a physical viewpoint this means that whilst small and large transverse momentum 
regions contribute equally to the BFKL equation, they are now weighted asymmetrically. 



The breaking of scale invariance in Eq. (|3.2|) allows one to regularize the collinear sin 



gularities in terms of e-poles. Once these poles have been factorized and subtracted, as 
in Eqs. (|2.36|) ,( p.38| ), one can perform the e — > limit and recover scale invariance at 



the expenses of an additional dependence on the factorization scale as dictated by the 
leading-twist behaviour of the parton densities (cross sections). 



A further consequence of the lack of scale invariance in the kernel of the gluon Green 
function is that no simple technique is available to diagonalize the integral equation (|3.2| ). 
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However, it is possible (see App. B) to solve it as a formal power series in as with 
^-dependent coefficients. The result reads as follows 



H 0) (k;/i,as,£)=5 (2+2e) (k) + 



r(i + e 



fvrk 2 



l+e 51 



fc=l 



£^(1) /k 2 \ £ 

N ~ £ T(l+~e) l~i^J 



(3.3) 



where we have explicitly introduced the MS-scheme factor S e = exp [— e + ln47r)], 

and the coefficients Ck are defined by the recurrence relation 



ci(e) = 1 , c fc+ i(e) = c*(e) 4(e) (fc > 1) 
with (r(z) is the Euler T-function) 



hie) 



i r 2 (i + £ ) 

e r(l + 2e) 



r(l + 2e)r(fce)r(l-fce) 
r((l + fc)e) r(l + (l-fc)e) 



-r(i + £ )r(i- £ ) 



(3.4) 



(3.5) 



The gluon Green function ( |3.3| ) is formally a distribution acting on the transverse- 
momentum space in n — 2 = 2 + 2 e dimensions. Following the procedure of factorization 
of collinear singularities described in Sect. 2.2, one can extract from the transition 
function T gg N as in Eq. (|2.36| ), thus leading to a 'renormalized' gluon density which, 
for finite e, is a very cumbersome /^-distribution. Therefore, it is more convenient to 
introduce the gluon Green function integrated up to the factorization scale /if, = Q 2 : 

G g °l N (a s (Q 2 /^r,e) = J d 2+2e k ^ 0) (k; a s , /i, e) 6(Q 2 - k 2 ) . (3.6) 

Note that does not depend on as e Q 2 /^ 2 independently, but only on the product 
asiQ 2 1 ' IJ 2 ) 6 ■ Performing the /^-integration of Eq. ( |3~3| ) we find 

oo 
k=l 



as q e 



T(l + e) 



Q 2 



ke 



(3.7) 



Introducing explicitly the parton indices, the collinear factorization in Eq. ( p. 21 ) reads 



G 



(0) 



ga, Af- 1 - ag,N- 



However, from the power counting in Sect. 2.3, we know that 



the quark anomalous dimensions 7^ n (and hence the corresponding transition functions 
r ga , n) are subleading at high energy. Therefore, to leading order in (as/iV) fc , we can write 

Gf g ] N (a s (Q 2 /f, 2 y,e) = G ggjN (a s (Q 2 /fM 2 y,e) T mN ia s iQ 2 /^) £ ,s) , (3.8) 



where ^ gg ,N is given by Eq. (|2.37 ), in terms of the gluon anomalous dimension 7 g9j 



N ■ 



By comparing Eqs. G3.7|) and ( [3.8|) , one can compute 7 gg ,Ar and the 'renormalized' Green 
function G 59 ,jv for any value of e. Moreover, from Eqs. ( |3.6| ) e Q2.36| ) it follows that the 
function Rn in Eq. Q2.38|) is related to the e — > limit of G g9t N 



Rwias) = G gg!N (a s (Q 2 /f/y,£)\ E=0 



(3.9) 



-16- 



The recurrence factor I n {e) behaves like 1/e for e — > 0. Therefore, in agreement with 
Eqs. Q2.37Q and Q3.8Q , the power series expansion (|3.7Q has at most a single e-pole for any 
power of as- More precisely, for e->0we have 



hie) 



ke 



and correspondingly 



Gfl N (a s ( Q 2 /n 2 Y,e) = 1 + £ 



~ exp 



las 

e N r(l + e) 



fc=i 

2\ £ 



1 a 5 e*^ 1 ) 

e77 £ r(i + £ ) 



0) 



(3.10) 



2\ £ 



M ( 1 + °< £: 



1 + 



OfifV 



(3-11) 



Comparing Eqs. ( |3.8|) e (|3.11[ ), we see that the slow departure of lfc(e) from its leading- 
pole approximation /^(e) — 1/fce gives rise to subleading collinear corrections, for both 
7 99j tv and i2jv, which are of relative order a| in the high-energy regime: 



799, N 



N 



O 



ivy 



iV 



i + c 



N J 



(3.12) 



Terms of relative order as and a| are present in the various Feynman diagrams con- 
tributing to the gluon Green function but they cancel in the sum. This cancellation, 
which is automatically embodied in the master equation (|3.2| ), does no longer occur in 
higher perturbative orders. The resummed expressions for 7 99i tv and Rn to the leading 
accuracy (as/N) n are derived in App. B. The results are the following. 

The dominant contribution 7^(05) to the gluon anomalous dimensions 

l gg ,N = lN{as) + 0{a s {a s /N) k ) (3.13) 
is obtained by solving the implicit equation (Fig. 6) 

l = |x(7»N) , (3.14) 
where the characteristic function xil) ls expressed in terms of the Euler ^-function 



X(7)= 2^(1) --0(7) --0(1-7) 



7 



l + £2C(2fc + l) 7 



2k+l 



k=l 



(3.15) 



and is the Riemann (^-function. The solution of Eq. (|3.14|) in power series of the 

coupling constant givesj] (C(3) ~ 1.202, ((5) ~ 1.037): 



71=1 



= 5>(£) =^ + 2C(3)(^) +2C(6)(^) +0((^ 



N 



N 



a s 



N 



a s 



N 



N 



(3.16) 



"The explicit values of the coefficients g n , up to n = 14, can be found in Ref. 
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The function R^(as) is given by 

( T(l - 7jv) x(1n) 



,r(l + 7jv) [-In x'(1n)] 



expW)+f V' (1) ^ (1 - 7) } , (3.17) 



where x an d x' are the characteristic function in Eq. (|3.15| ) and its first derivative, respec- 
tively. The ag-dependence of the r.h.s. in Eq. ( |3.17| ) is implicit in that of the gluon anoma- 
lous dimension 77V = 7jv(as). The first perturbative terms are (C(4) = 2£ 2 (2)/5 ~ 1.082): 

^K)- + ?C(3,(f) 3 -^,(|) 4 + fc(5,(f) 5 + ((f) 6 ) . ,3,S) 
Some comments are in order. 



The result in Eqs. ( |3.13| ) , ( |3. 14j ) for the gluon anomalous dimension is exactly the cele- 



brated BFKL anomalous dimension [T^| . In the present paper this result has been derived 
by consistently carrying out the procedure of factorization of the collinear singularities in 
dimensional regularization. However it is worth noting that the same expression for the 
anomalous dimension can be obtained by using alternative and less sophisticated regular- 
ization prescriptions of the collinear singularities (for instance, considering the Eq. ( |3.2|) 
directly in n = 4 dimensions and regularizing it by keeping the incoming gluon slightly 
off-shell, via the replacement 5^(k) — > 5(k 2 — /i 2 )/vr in the inhomogeneous term p4]j35[| ). 
The reason for this has to be traced back to the property of the kernel in the master equa- 
tion ( |3.2| ) of being collinear regular for e — > 0: the collinear divergences in the solution 
of Eq. (|3.2|) originate only from the fact that the inhomogeneous term is not sufficiently 
smooth for k — > 0. As a result, in the small- N limit, the gluon anomalous dimensions 
to leading accuracy (as/N) k are regularization/ factorization scheme independent (and 
related via Eq. fl3.14| ) to the eigenvalues xil) °f the BFKL kernel) within a wide class of 
schemes. This class includes all the schemes which do not introduce pathologically sin- 
gular terms of the type a k s /N k+p (p > 1) in the perturbative calculation at high energy 
[], or, more precisely, which do not violate the high-energy power counting discussed in 
Sect. 2.3. 

The BFKL anomalous dimension in Eqs. ( |3.14|) , (|3.16|) departs rather slowly from its 
one- loop contribution. However for very small values of x, corresponding to sizeable 
values of a s /N ~ 0(1), 7^ increases quite fast and for a s /N = (4 hi2) _1 reaches the 
saturation value 7 = 1/2 at which x(l) has a minimum (Fig. 6). For still larger values of 
as/N there are two complex conjugate branches of 77V, coming from the pinching with 
the symmetrical solution of Eq. ( |3.14|) at 7 = 1 — 7^. Therefore the resummation of the 



singular terms (as/N) k builds up a stronger singularity at iV = A = 4as In 2 ~ 2.65 as 
|i"2| . As discussed in Sect. 1, this branch point singularity (known as the perturbative 



QCD pomeron) is responsible for the steep behaviour x A of the gluon density f g (x, /i 2 ) 
generated by the perturbative QCD evolution at high scale /i 2 . 

^This feature of the gluon anomalous dimensions was first pointed out by T. Jaroszewicz |36f . 
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The function Rw(as) in Eq. ( |3.17| ), on the other hand, depends on the factorization 



scheme more than the gluon anomalous dimensions. For instance, going from the MS- 
scheme result ( |3.17 ) to the MS scheme (in this scheme the transition function T gg is 



obtained by setting S e — 1 on the r.h.s. of Eq. ( |2.37| )), RN{as) has to be multiplied by 
the factor exp [— 7^(0:5) (ip{l) +ln47r)]. This scheme dependence of Rn(ccs) has to be 
compensated in physical observables by subleading contributions of order a>s(o>s/N) k in 
the anomalous dimensions. 

Since R^{as) is related to the e — > limit of the renormalized gluon Green function 
G gg ,N (see Eq. fl3Jp ), it can be regarded as the normalization factor of the perturbative 
QCD pomeron. Note, in particular, that Rjsr(as) is singular at the saturation value 
7 = 1/2 of the BFKL anomalous dimension 7^(05): 

R N (a s ) ^ const. (- 1 j , ( lN -+ 1/2) . (3.19) 

This behaviour is related to the branch point singularity of the BFKL anomalous dimen- 
sion at N = 4 a>s In 2 and signals the ultimate failure of the leading-twist approach to the 
mass singularity factorization and the onset of the multi-Regge factorization at extreme 
energies (see Ref. || for a more detailed discussion of this issue). 

We have so far considered the leading high-energy behaviour of the gluon Green func- 
tion Qfg- In order to complete our discussion on the gluon channel, we conclude this 
section by examining the Green function Qj® with an incoming quark or antiquark. 

Let us denote by the fc^-dependent Green function defined by replacing Qj® by 
Qj® on the r.h.s. of Eq. ( p.35| ). The corresponding master equation to leading accuracy 



(as/N) k is the following 

FlUk.;as,n,e) = 77-7:77- — 777 ~~ri + 



" - vV W ' " C A N (27r/i) 2£ vrk 2 N J (27T/i) 2e vrq 2 

(°) (v _ ~. n ,, *\ _ k ■ ( k - q) tt(o) 

(k-q) 



•^^(k-q; a s,V,e) - -77- vj- -^°M k ; atf.^e) L(3.20) 



Comparing Eq , (|3.20p and (^.2[), we see that T gN an d •7 r /v ; fulfil the same integral 
equation, apart from a different inhomogeneous term. In particular, it is straightforward 
to check that the solution of Eq. (|3.2(J|) can be expressed in terms of the pure-gluon Green 



-(0) 



function JF^ as follows 



a s , H,e) = %Z M\k; a s , fi, e) - 5^\k) 



C 



The fcj_-integrated Green function G g °J is thus given by 

G g °l N (a s (Q 2 /n 2 y,e) = jd 2+2£ k ^%(k; a s , /*, e) 6(Q 2 - k 2 



Cf 
C a 



G 9 ° g ] N (a s (Q 2 /^y,e)-l 



(3.21) 



(3.22) 
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The result (|3.22|) allows us to compute the anomalous dimensions r y gq ,N by per- 
forming simple algebraic manipulations. We have first to factorize from N the 
collinear singularities according to Eq. ( j2.21| ). Introducing explicitly the parton indices, 

Ggq,N = J2 a G g a,N^aq,N, and using the fact that the quark transition functions T q b are 
subleading at high energy, we obtain: 

G { 2 N (a s (Q 2 /^y,e) = G gg , N (a s (Q 2 /^r,e)T gq , N (as(Q 2 /^y } e) 

+ G gqtN (a s (Q 2 /fi 2 ) £ ,e) . (3.23) 

Inserting Eqs. fl3~23l) and Q) into Eq. ( ^22]) we get the identity 

Ga [Gg gi N(®si £j r 9gi iv(a5, e) + G gqiN (a s , £)] = C F [G ggjN (a s , e) r g9iA r(a5, e) — 1] , 

(3.24) 

which we can rewrite in the following form: 

[Ca G gq .N{&s, e) — G F (Ggg 7 ^(as, e) — 1)] — - 

= C F [T ggtN (a s ,e) - 1] - C A T gqiN (a s ,e) . (3.25) 

We now notice that, order by order in perturbation theory, the renormalized Green func- 
tions G a b are regular for e —>■ whilst the transition functions Tab are series of poles. 
Therefore the only solution to Eq. fl3.25|) is 

C F 

G gq , N {ois,e) = — [G ggtN (a s ,e) - 1] , (3.26) 
C 

r 5<? ,Ar(a 5 ,e) = — - [T ggtN (a s ,e) - 1] . (3.27) 

^A 

In particular, since to leading accuracy (as/N) k we have 

1 rots s e dot 

F gq ,N(as,z) = ~ J Q — hgq,N(a) + 'jgg,N(a)T gqjN (a,e)} , (3.28) 

1 f a s da 

Fgg,N(as,z) =1 + ~ — lgg,N(a)T gg!N (a,e) , (3.29) 

& J Gi 



from Eq. (|3.27 ) we obtain: 



C 

7 9q ,N(as) = 7 ^7gg,N(as) + 0(a s (a s /N) k ) 



C A 

( F i N {as) + 0(a s {a s /N) k ) , (3.30) 



C A 

1n{ols) being the BFKL anomalous dimension. We see that the coefficients of the leading 
terms (as/N) k in 7 9a ,iv (a = Q, g) are equal, apart from an overall factor given by the 
ratio of the colour charges of the initial-state parton a. 
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4. The quark Green functions 

From the high energy power counting discussed in Sect. 2.3, we know that the gluon 
channel dominates to leading logarithmic accuracy in the small-x regime. Beyond lead- 
ing order, however, the quark sector has to be considered on an equal footing with the 
gluon sector. Actually, from a phenomenological viewpoint, the knowledge of the next-to- 
leading quark anomalous dimensions may be more relevant than that of the corresponding 
corrections to the gluon anomalous dimensions. The reason for this is that the most ac- 
curate information on small-x parton densities is coming out from HERA data |J7| on 
deep inelastic structure functions, which couple directly to quarks (and not gluons). 

In order to compute the (flavour-singlet) quark anomalous dimensions we consider the 
Green function Qfj$ : 

G { 3 = Y, K{ $Gto! > (4.i) 
b 

or, more precisely, its expansion in 2GI kernels (Fig. 4b), and we apply the high-energy 



factorization formula ( 2.29 ). We thus arrive at the analogue of the fc^-factorization 
formula (gjj : 



C" W)= l<?~ 2 *l - [^J a %,k)^pi ^\y,k-as^,e) , (4.2) 

K / | k=yp+k± 




where (Fig. 7) is the 2GI kernel Kf® to the lowest order in as (//, v and a, (3 

are respectively the spin indices of the incoming gluon and outgoing quark). Moreover, 
since we are interested in the anomalous dimensions rather than the coefficient function, 
it is also convenient to apply the collinear projector Vc and consider (as in the gluon 
sector) the Green function G q °J integrated over and the invariant mass q 2 up to the 
factorization scale \j? F = Q 2 (we use the Sudakov parametrization q^ = x p^ + q^_ + (q 2 + 
q 2 ) p' 1 /2xp ■ p ) 



dq 2 d 2+2e q ( fi \ 



G$ aP (q,p) e Q 2 -|g 2 |) . (4.3) 



2(2tt)*+*« \2p.q) a /™ 
Inserting Eq. ( fOj) into Eq. (|4.3|) and taking the iV-moments we get 
Gfl N (a s (Q 2 /» 2 y,e) = J d 2 + 2 ^ K qgtN (k 2 /Q 2 , a s {Q 2 /^Y;e) ^(k; a s , //, e) . (4.4) 



The factorization formula ( |4.4j ) relates the quark Green function to the /^-distribution 
J 7 ^ in the same way as Eq. fl2.34j ) relates the heavy-flavour cross section to the gluon 



density. The off-shell kernel K q9t at, obtained from K^ \q,k) after integration over q and 
explicitly calculated in App. C, is indeed regular for iV — > so that the resummation of 
the singular terms (as/N) h in G q °g N is achieved by fc^-integration of the corresponding 
terms in J 7 ^. Note, however, an essential physical difference with respect to Eq. (|2.34j ). 
Unlike the hard cross section &n, the off-shell kernel K qg ,N is not collinear safe. The 
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collinear singularity arises from the integration of the splitting process g — > qq in Fig. 7b 
and shows up as an e-pole in the n-dimensional on-shell case k 2 = : 

■ -Jn , (4.5) 

or, alternatively, as a logarithmic divergence in the on-shell limit k 2 — ► for the four 
dimensional case: 



K qg (z,k 2 /Q 2 ,a s ,e = o) = z||e(Q 2 -zk 2 ) (4.6) 

((l-z) 2 + z 2 ) ln^y-(l-6z(l-*)) (l 



2tt 
Tr 



zk 2 ^ 



Q 2 



Note that the coefficients of both the e-pole and ln(Q 2 /k 2 )-term are correctly proportional 
to the n-dimensional Altarelli-Parisi splitting function Pj® in one-loop order: 

P^s ) = T ^-^f + £ . (4.7) 

Actually, the full off-shell kernel K qg ^ in Eq. (|4.4|) can consistently be interpreted as 
the integral of a generalized off-shell (and positive definite!) splitting function Pig ■ The 
explicit calculation in App. C gives 



k qg {z 1 k 2 /Q 2 ) a s {Q 2 /^ 2 y ] e) = zQ(Q 2 -zk 2 ) [ 

J 



(l-z)(0»-*k») dq 2 /q2\ £ 

q 2 U 2 , 



where the splitting function is 
P£\z,k 2 /q 2 ;e)=T R (— 



q 2 



f^T^^kVqV), (4.8) 



'l-z) 2 + z 2 + e , , 2n , 2 k 2 ^ 



q 2 + z (1 - z) k 2 



+ Az z (l-z 



1 + e ' q 2 



(4-9) 

and q = q— z k denotes the boost-invariant (along the /c-direction) transverse momentum 
transferred in the splitting process g — > qq of Fig. 7b. 

The fact that the kernel K qg is not collinear safe prevents us from taking the e — > 
limit in Eq. (|4.4| ) even after having factorized the transition function T gg ^N as in 
Eqs. ( |2.34|) , (|2.36|) . The non-polynomial e-dependence of K qg (due to both P« and the 
phase space integration in Eq. ( ^.8|) ), as much as the off-shell dependence of P qg \ are 
thus responsible for non-trivial (and factorization scheme dependent) quark anomalous 
dimensions 7 gg ,Ar in higher perturbative orders. 

Using the power series solution ( |3.3| ) for the gluon distribution J-"ff and performing the 
/^-integration in Eq. ( f4.4|) , we can determine the corresponding power series expansion for 
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the quark Green function G q °J jN . Considering the N — > limit for K q9l N (i-e., neglecting 
subdominant corrections at high energy) we find (see App. C) 



'Q 2 



4 + e 



2tt 



r(l + e) \fi 2 J e (3 + e)(2 + e) 




as 9 e^g) 
AT e r(l + e) 



^dfc(e) 1 + O (<4(a s /iV) fc ) , (4.10) 



where the coefficients dfc(e) are expressed in terms of Ck(s) in Eq. (|3.4|) as follows 
1 4 + (1 - 3k)e T(l + ke) T(l -ke) T(l + e) T(4 + e 



d h (e) 



k + 1 



r(l + (l + fc)e)r(4+(l-fc)e) 



(4.11) 



On the other side, the procedure of factorization of the collinear singularities in 
Eq. flUO) gives Gqg jy — Gqa,N^ag,N> which, neglecting terms of order a^^as / TV)' 5 , 
can be written as 



qg,N 



= G q9;N (a s (Q /// ) £ ,e) T ggiN (a s (Q j \x 
+ r qgiN (a s (Q 2 /n 2 ) £ ,e) , 

1 f a s Se da 



2\e 



(4.12) 
(4.13) 



Remember that G qg ^N is finite for e — > order by order in as- The structure of the e-poles 
on the r.h.s. of Eq. (|4. 12j ) thus defines G q9t N and the transition function T q9t N uniquely. 
It follows that, comparing the expansion ( 4.10|) with Eqs. (|4.12|) ,( |4.13l ) and using the 
known result for T g9jN (i.e. the leading-order gluon anomalous dimension ( |3.13|) ), we can 
compute the anomalous dimension 7 g9i Ar (and the renormalized Green function G qg ^) 
order by order in perturbation theory. 

We have explicitly performed this calculation up to the sixth order and the result for 
the quark anomalous dimensions in the MS scheme reads 



7qg,N(as) 



2tt R 3 } 



5 05 14 ras_Y 
3 N 9 \NJ 



+ 



g + 2C(3) 



NJ 



+ 



25 
~6~ 



C(3) 



NJ 



+ 



146 14 

^ + T C(3) + 2 C (5) 



122 
243 

as 
N 



i R — 

2tt 3 



1 + 1.67— + 1.56 

N 



atsY 
NJ 



N 



+ 



as 
N 



(4.14) 



The calculation is straightforward but already very cumbersome to this order in pertur- 
bation theoryf], and we are not able to provide an explicit resummed formula for 7 g9 to 



*Note that, because of the cancellations involv ed in the gluon sector (see Eqs. ( 3.10 ),( p2| )), the first 
three coefficients in the perturbative expansion ( 4.14 ) are entirely related to the e-dependence of the 
off-shell kernel K qg , n ■ 
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all orders (some all-order features of 7 gg are presented in App. C). However we emphasize 
that this is just an (open) algebraic problem related to the use of the MS factorization 
scheme, that is the highly non-trivial e-dependence of the coefficients dk{s) in Eq. ([4.1 1|) 
and the non-local structure of the e-poles in Eq. ( |4.13|) . The series ( fj.lOj ) for the quark 
Green function indeed contains all the necessary information on the anomalous dimen- 
sions 7 gs ,v to any perturbative order a>s(cts/N) k . In the next section we show that, 
choosing a different factorization scheme of mass singularities, we can explicitly resum all 
the next-to-leading terms a>s(o>s/N) k in 7 g9 ,v- 



The coefficients of the first two terms in the curly bracket of Eq. Q4.14|) agree with 
the known one- and two-loop anomalous dimensions in the MS scheme pTp8| . The 
higher-order contributions can be used to estimate the effect of these small- X corrections 
at intermediate values of x. In particular, the 0{a%) term in Eq. ( |4.14| ) can be combined 
with the existing 0(a 2 s ) calculations of the coefficient functions for the DIS and 
Drell-Yan |JT[ processes in order to check the stability of the fixed-order perturbative 
expansion in the x-range accessible at present. Note also that, unlike the case of the 
gluon anomalous dimensions ( |3.13| ),( |3".16| ), all the perturbative coefficients in Eq. ( |4.14| ) 
are non- vanishing (and positive). Therefore in the quark sector one may expect |22| (and 
finds [[H|) an earlier departure from the fixed-order perturbative behaviour. 

The flavour- singlet anomalous dimension 7^ N (App. A) starts in order a 2 s and its 
resummed expression at high energy is related in a simple way to that for 7 gff ,7v, namely 



lqq,Ni a s) 



C a 



i qg M a s) - -^t r - 



+ 0(a 2 s 



(a s /N)> 



(4.15) 



The derivation of Eq. ( [4.15[ ) is similar to that of the analogous Eq. Q3.30| ) for the gluon 
sector. One has to apply the high-energy factorization procedure of Sect. 2.3 to the Green 
function G qq ^ (here, the superscript 5* denotes the fact that the incoming and outgoing 
quarks carry a different flavour), thus obtaining the following fc^-factorization formula 



G s q ^(a s (Q 2 /^ 2 r,e) = J d 2+2s k K mN (k 2 /Q 2 , a s {Q 2 / ^f-e) J*%(k; a a ,n,e) . 

(4.16) 

This equation differs from Eq. (|4.4j) only by the replacement of the pure-gluon k±- 
distribution by J 7 - N (see Eq. (|3.20|) ). Therefore, using Eq. ( |3.21| ), one can relate 

the two Green functions G qq ^ and G q °J : 



(o) / ( Q 2 

G qq,N \*S [y 



C a 



e)-K 



19, N 



0,a s 



'91 

y 



17) 

The factorization of collinear singularities in G qq °x is achieved by the following high- 
energy relations 



r s (0) 
qq,N 



qq,N 



^qq,N( a s, + G qq N (as, e) + v(«5, £) r OT ,v(a5, s) 



(a s ,e) 



1 

£ JO 



a 



lqq,N{ a ) + lqg,N{a) ^^(a, s) 



(4.18) 
(4.19) 
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Using Eqs. ( |4. 18| ) , ( f£T§| ) , ( |3.27| ) one can express G^ -* as a function of G qql 7^, G qg , 7 q9 , 



T gg . Analogously, by means of Eqs. ( }4.12[ ), (|4.13| ) one can express G q °J as a function of 
G qg , 7 99 , T gg . Inserting the expressions obtained in this way into Eq. (|4.17| ), one ends up 
with the equation 

r a sS E da 



G qq,N( a S,£) + 



e Jo 



■7, 



5 

qq,N 



at 



C a 



G 



qg,N(0ts, e) + ~ 



a: 

1 r°*sS e da 







— lqg,N{pt) - K q g :N (k = 0, a s ; e) 



a 



(4.20) 



from which the result ( 4.15|) follows. 



It is worth noting that the all-order relation (|4.15|) between 7 ( 5 



as its leading-order analogue (|3.30 



qqN and j gg ,N, as well 
have been derived using only algebraic identities 



at high energy, with no reference to the details of the dimensional-regularization pre- 
scription. In particular they remain true using dimensional reduction 



a dimensional 

regularization scheme which explicitly respects supersymmetric Ward identities. In this 
sense, Eqs. ( |3.3U| ) an d ( 4.15 ) can be considered as high-energy limits of the supersym- 
metry identity 7,^ + 7 OT = 7 g9 + ^gg valid in iV = 1 supersymmetric Yang-Mills theory, 
i.e. for Cp — Tr = Ca- It follows that in the supersymmetric case the gluon anomalous 
dimensions 7 gg and 7 OT coincide also at the next-to-leading level as(as/N) k . This prop- 
erty may be useful as a technical tool to check and simplify the calculation of the (still 
unknown) next-to-leading corrections in the gluon sector. 



5. Deep inelastic scattering at small x 



5. 1 Structure functions and parton densities 

The cross section for deep inelastic lepton-hadron scattering is given in terms of the 
customary structure functions Fi(x,Q 2 ) (i = 1,2,3). Here Q 2 denotes the square of 
the momentum transferred by the scattered lepton and x is the Bjorken variable. In 
the following we only consider the scattering process occurring through the exchange of 
a single photon. As far as the hadronic component is concerned, this approximation 
simply amounts to neglecting F 3 , which is a non-singlet structure function and, hence, 
non-singular at small x (see Sect. 2.3). We also present our results in terms of the struc- 
ture functions F 2 and Fl, Fl(x,Q 2 ) = F 2 (x,Q 2 ) — 2xFi(x,Q 2 ) being the longitudinal 
structure function. 

In the naive parton model the DIS structure functions are related to the parton den- 
sities of the incoming hadron as follows (e^ are the quark charges) 

Nf 

F 2 (x,Q 2 )=J2e 2 {fqSx,Q 2 ) + k(^Q 2 )} , (5-1) 



i=i 



F L (x,Q 2 ) = . (5.2) 
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Taking into account perturbative QCD corrections in leading-twist order and using the 
factorization theorem of mass singularities, Eqs. (|5.1|) , flO|) become 




Ffa Q 2 ) = i- I V e\ I if (x, Q 2 ) + i^ 5 (x, Q 2 ) , (i = 2,L) , (5.3) 



where the singlet and non-singlet components F s and F NS are given by 
Ff(x,Q 2 ) = I' — [cS{z-a s ^ 2 F ),Q 2 /ti)~f s {x/z^ 2 F ) 

J x Z 

+ Cf(z;a s (fi 2 F ),Q 2 /fi 2 F ) f g (x/z^ 2 F )] , (5.4) 



?NSi 



(x,Q 2 )= f —C^ s {z ] a s ^ 2 F \Q 2 ^ 2 F )Y.e 2 f^\x/z^ 2 F ) . (5.5) 

Z 3=1 

Here fs and are respectively the singlet and non-singlet quark densities defined 

in App. A, and Cf (A = S, g, NS) are the coefficient functions^ computable as power 
series in as {up)- Note that, because of the naive parton model relation ( |5.1|) , we have 
C£{z) = 6{l-z)+0{a s ) for A — S, NS, and C 9 2 {z) = 0{a s ). Therefore also the gluon 
density f g contributes to F 2 beyond the lowest order. Moreover, since C F = 0(as) is 
not vanishing, the Callan-Gross relation (|5.2|) is violated. 



At present, the coefficient functions Cf are completely known up to 0(a s ) |p0,p9|. 
The non-singlet coefficients Cf s are not enhanced by In a; terms in higher orders. The 
all-order resummation of the logarithmic contributions a^ +2 ln n x (or as(as/N) n+1 in 
the A^-moment space) to Cf and Cf is performed in the next subsection. 



5.2 Coefficient functions 



The coefficients functions Cf are evaluated starting from the expression Fi = 
52a Fia fa°^ °f the hadronic structure functions Fj in terms of the partonic structure 
functions F-® , and then performing the collinear factorization as in Eq. fl2.15|) . Consid- 
ering first the gluon structure functions F^\ we have 



7 (o) 



1 



N, 



ir. £ 



N 



CfT gg + Cf2N f T 



f L 19 



(5.6) 



On the other side, from the high-energy factorization in Sect. 2.3 (see Fig. 4a), we obtain 



2\£ 



/ \i=i 

"i dz 




2+2ei 



(5.7) 



f — a!(z,k 2 /Q 2 ,a s (Q 2 /fi 2 Y;e) F$\x/z, k; a s , //, e) 

J x Z 



*The coefficient func tions in Eqs. (5.4),( |5.5| ) are normalized according to our notation in Sect. 2.1 (see, 
in particular, Eq. ( 2.11 )). This normalization differs from that often used in the literature [ j30| . 
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where J-^ ' is the fcj_-dependent gluon Green function in Eq. ( |2.35| ), and of are obtained 
by applying the high-energy projector Vh to the lowest-order qq contribution to the 
1*9 1*9 absorptive part A^ u (Fig. 5), as in Eq. ( |2.30|) . The off-shell cross sections of 
have been explicitly evaluated in Refs. |8j (i = 2) and [32| (i = L) for the case of massive 
quarks and n = 4 dimensions. The generalization to the massless case and n = 4 + 2 e 
dimensions is straightforward. 

Eq. ( |5.7| ) is the dimensionally-regularized version of the fc^-factorization formula ( |1.2| ) . 
We can thus properly address the issue of collinear-singularity factorization. The main 
point to be noticed is that the off-shell cross sections a 9 L and of have a different collinear 
behaviour. As a consequence of the Callan-Gross relation ( |5.2| ) (i.e., the fact that on-shell 
partons do not couple directly to a longitudinally polarized photon), a g L is collinear safe. 
Its on-shell and e = limit is indeed finite: 



a 9 L (z,k = 0,as(Q 2 /fi 2 ) £ ;e) = -^S £ ^ N f T R 8 z 2 (1 - z) + 0(e) . (5 



On the contrary, af is not collinear safe and its on-shell limit has an e-pole proportional 
to the lowest-order Altarelli-Parisi splitting function P^> in Eq. Q4.7|) : 



&l(z, k = 0, a s (Q 2 /^Y; e) = ^ S £ (^) N, \ [zPjg>(z; e) + 0(e) 



(5.9) 



Therefore, the factorization structures ( |5.7|) for F^g and F^J are respectively analogous 
to Eq. (|2.34|) for the heavy-flavour case and Eq.( [4.4|) for the quark Green function. The 
factorization of collinear singularities has to be carried out accordingly. 

We start considering the longitudinal structure function. In the high-energy limit, 
both the coefficient function Cf and the transition function T qg are of order ots{ois/N) k 
{k > 0). Therefore from Eq. (|5.6p we have 



p(o) 




Clr 9g + 0( a 2 s (us/N) k )} , (5.10) 



f 

and we see that, in order to compute Cf , we have to factorize T gg on the r.h.s. of 
Eq. ( |5.7|) . Using the factorization formula ( |2.36 ) and proceeding as in Sect. 2.4, we arrive 



at the following expression for the iV-moments of the longitudinal coefficient function 

C1tv(«5,Q7/4) = Kn (lN(a s )) R N (a s ) (Q 2 / '/4) w(as) + 0(a 2 s (a s /N) k ) , (5.11) 
where the function Iil,n{i) is given in terms of the off-shell cross section <j£ : 

h L ,N(l)=lJ o Kn(^/Q 2 ^s;b = 0). (5.12) 



The explicit evaluation of the function h LtN (^) requires the knowledge of the off-shell 
cross section a 9 L in n = 4 dimensions. The analogous calculation for massive quarks has 
been performed in great detail in Refs. |8|j32]|. Therefore in this paper we limit ourselves 
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to presenting only the final result. In particular, since the off-shell cross section a 9 L 
in Eq. ( |5.7|) vanishes uniformly in k± in the high-energy limit £ — *> (as discussed in 
Sect. 2.3, this property follows from the fact that only 2GI kernels contribute to a 9 L ), the 
function h F ,N(l) is weakly iV-dependent: 



h L Al)= h L { 7 )(l + 0{N)) 



N 



(5.13) 



In order to evaluate the dominant terms as (as/N) k in Eq. ( |5.11| ), we can thus set N = 
in Eqs. ( |5.12|) , (|5.13|) . Computing explicitly the corresponding function h F {pf) , we find 



h L ( 7 ) = ^N f r R 4(1 ~ 7) r 3 (i-7)r 3 (i + 7 ) 

ULU) 2ti f R 3-2 7 T(2 - 2 7 ) T(2 + 2 7 ) ' 



(5.14) 



The results in Eqs. ( |5.1l] ),( |5.14|) give the resummed expression for the coefficient func- 
tion CI N in the MS scheme to the logarithmic accuracy as (as/N) k , including the 
corresponding dependence on the factorization scale /j,p . The resummation effect is in- 
corporated in ( |5.11| ) through the (as/iV)-dependence of the BFKL anomalous dimension 
(|3.14|) and the 7-dependence of and h L N as given by Eqs. (|3.17| ) and (|5.12|) . 



Let us now consider the structure function F 2 . As noted above, the pattern of collinear 
singularities in Eq. (|5.7|) is similar to that in the corresponding Eq. (|4.4j ) for the quark 
Green function. Eq. (|5.7| ) thus contains all the relevant information on both the DIS 
coefficient function C 2 and the quark anomalous dimensions 7 gg . As a matter of fact, 
since in the small- N limit C 2 N = 1 + 0(as(as/N) k ), from Eq. ( |5.6| ) we obtain the 
analogue of Eq. (|4.12j ) : 



F 2 {0 U»s(Q 2 /^Y,e) 



2\£ 




Cl N {a s {fi F /fi ) e ,Q /fj. F ;£)T gg>N (a s {fJ, F /fj, ) £ ,e) 



+ 2N f r WlN {a s (v F /ii 2 ) e ,e) + 0{a 2 s {as/N) 



(5.15) 



Starting from Eq. ( |5.7|) and factorizing the collinear singularities according to Eq. (|5.15| ), 
one can compute Cf N and 7 g9i Ar order by order in perturbation theory. 

Obviously, this algebraic problem is by no means simpler than that encountered in 
Sect. 4 for the evaluation of 7(?g ,Ar and we are not able to provide an explicit resummed 
formula for C 2 n- Nonetheless, Eq. (|5.7|) can be used in a very simple way for deriving 



an all-order relation between C 2 N and l qg ,N- 



The main observation is that, setting fi 2 F = Q 2 on the r.h.s. of Eq. ( |5.15| ) and 
performing the derivative with respect to Q 2 , we obtain a factorized structure similar to 

Eq. (prig) : 







2\e 



l 99 ,N(a s (Q 2 /fi 2 y)Cl N (a s (Q 2 /v 
d 




(5.16) 



+ ea s 



da s 



a N (a s (Q 2 /v 2 y,l;e) 



l;e) + 2N flqg , N (as(Q 2 /fi 2 ) £ ,e) 
T gg , N (a s (Q 2 /li 2 y,e) + 0(a 2 s (a s /N) h ) 
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Here, all the collinear singularities are factorized into the gluon transition function T gg 
whereas the term in the square bracket, much like C 9 L in Eq. ( |5.10| ), is finite as e — > 0. 
Stated differently, the off-shell kernel da 2 /d\nQ 2 in Eq. ( |5.7|) (although not erf itself) 
is collinear safe, and its e = limit is related to the linear combination 7 5g C| + 2Nf / ~f qg . 
Therefore, proceeding as in the case of the longitudinal structure function, we obtain 

77v(a s ) C$ N (a s , Q 2 /^ 2 f = 1) + 2N f l qg ,N(a s ) = h 2jN {jn{(Xs)) Rn((*s) 

+ 0(a 2 s (a s /N) k ) , (5.17) 

where is given in Eq. ( |3.17| ), 7^(0:5) is the BFKL anomalous dimension and the 
function /2.2,7V (7) * s 

r°° dk 2 ( k 2 \ 7 d 
h,Al)=lj ^(^) ^^(kVOW-O). (5.18) 



The off-shell kernel dd\ N /d\aQ 2 in n = 4 dimensions has been computed in Ref. 



for the case of massive quarks. Performing the massless limit, we find 

h 2 , N h)= h 2 (l)(l + 0(N)) , (JV->0) , 

3 7 2 ) r 3 (i- 7 )r 3 (i + 7 ) 
r(2-2 7 )r(2 + 2 7 ) 



h( \ aS AT T 2 ( 2 + 3 7 

h 2 { 1 ) = -N f T R — 



(5.19) 
(5.20) 



The results in Eqs. (|5.17j ),( |5Tl9"l) ,( |5.2d|) provide the explicit resummation of the contri- 
butions ots(ots/N) k to the coefficient function C 2N in terms of ^(tjv) an d the quark 



anomalous dimensions 7 gg ,jv. In particular, using the six- loop expression ( |4.14|) for 7 gs ,Ar, 
it( 

ization scale \j? F is given by 



one can compute C 2 N up to the five-loop order. The dependence of C 2 N on the factor- 



C 2 M(a s ,Q 2 /f4 



+ 2N 



C% N {ocs,Q 2 lv F 

lqg,N{0Cs) 



7iv(os) 



1n(ocs) 



Q 



2 \ 7jv("s) 



+ 0{a 2 {a s /Nf) .(5.21) 



The singlet coefficient functions Cf , C 2 in Eq. (|5.4j) can be evaluated starting from 
parton structure functions F^ g with an incoming quark. These structure functions fulfil 
a fc^-factorization formula similar to Eq. (|5.7|) with the replacement — > 
being the fcj_-distribution with an incoming quark in Eq. ( |3.20D . One can apply the 
algebraic manipulations analogous to those used in Sects. 3 and 4 for evaluating 7ot and 



7 £ 9 , thus obtaining the following colour charge relations 



r<ps 
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+ O (a 2 (a s /N) k ) 



(5.22) 



2 1n^) 



+ (a|(a 5 /iV) fc 

(5.23) 
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Here, for the sake of convenience, we have introduced the pure-singlet coefficient functions 
Cf s = Cf — Cf s . They have the same singular small-x behaviour as the singlet functions 
Cf and start in 0(a 2 s ) in perturbation theory. 



The perturbative expansions of the resummed results derived in this section read: 

2 r 40 



Cl N {a Sl Q 2 //4 = 1) = ^T R N f Ul- ^ + 



as 
N 



+ -C(2) + -C(3) 



+ o 



N J 



+ 



^-fc( 2 )^C(3)- 6 C( 4 ) 




^T R N } \ { 1 - 0.33| + 2.13 (|) 2 + 2.27 (f )' + 0.43 



+ o 




N 



N J 



27 



(5.24) 



Cl N (a s , Q 2 //4 = 1) = ^-T R N f \\ 1 + 



f-2 C (2) 



1234 -H<(2) 



51 




Z«?_n ((2) + » ((S )-l2C(4) 
243 9 SV 7 9 SV 7 SV . 



aJ 



910 

C(2) + ^C(3) 



28 24 
y C(2)C(3)-26C(4) + y C(5) 



(5.25) 



2 -TnN } \\l + lA^ 4 0.71 



16.43 r^v 



16. 13 (^j +39-11 (I)' 




The first two coefficients in Eqs. ( |5.24| ),( |5.25| ) agree with those recently computed in 
Refs. [30, l0||. We regard this agreement as a non-trivial check of our results. Note also 
that the three- and four- loop coefficients (and the five-loop coefficient in Eq. (|5.25|) !) 



are substantially larger than the two-loop ones. We thus argue that the higher-order 
contributions computed in this paper may have a phenomenological relevance already at 
the values of x accessible at the HERA ep-collider 



Concluding this subsection, we point out that /c^-factorization formulae similar to 



Eq. (5.7) have recently been used [19| with the phenomenological aim of relating the 
original BFKL equation |12| to the DIS structure functions. 



5.3 The DIS factorization scheme 

In the previous sections we have repeatedly noted that the parton densities are not 
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physical observables. Indeed they depend on the regularization/f actor ization scheme used 
for removing the parton level collinear singularities. This freedom in defining the parton 
densities means that, starting from the MS densities /„, one can introduce a new set f' a 
of parton densities via an invertible transformation 

fa,N(v 2 ) = 5> a6 ,ivMM 2 )) iU(/i 2 ) • (5.26) 
b 

Obviously a similar transformation applies to the coefficient functions, in order to leave 
the physical cross section unchanged. The evolution of the new parton densities with /j, 2 
is controlled by the new anomalous dimension matrix 

(5.27) 

ab,N 



lab,N 



d 



T-l 



The transformation matrix U a b has a power series expansion in as such that U a t, t n{&s) 
5 a b + 0(as), and has to fulfil the following physical constraints : 

i) flavour and charge conjugation invariance 

Uq t g = Ug ig = U qg , U gqi = U g q i = U gq 

U mj = U m = - i-) U™ + JL (U^ + U SS ) , (5.28) 



Uqtqj — Uq iqj 



ii) fermion number conservation 

Uq iqj) N=Q — U qi g Jt N=0 = ^ij (5.29) 

or, equivalently, U»? N=0 = U»? N=0 = l/j&fa 

iii) longitudinal momentum conservation 

E Uab,N=l = l ■ (5-30) 



a 



Eq. fl5.28| ) is the analogue of Eq. ( |A.4| ) in App. A for the anomalous dimensions 7^. In 
particular, it guarantees that the flavour singlet and non-singlet sectors are decoupled in 
any regularization/f actorization scheme. The matrix components U^ s , , U^ v ' intro- 
duced in Eq. ( |5.28|) act on the flavour non-singlet parton densities, whilst Uss, U qg , U gq , U gg 
control the transformation on the singlet sector. 

Higher-order QCD calculations for hadron collisions are usually performed in two 
different factorization schemes of collinear singularities, the MS scheme, used so far in this 



paper, and the DIS scheme [41| . After having regularized the collinear singularities in the 
parton matrix elements, the DIS-scheme parton densities fa DIS ^ are defined by enforcing 
the constraint that the DIS structure function ^(x, Q 2 ) has the same expression as in 
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the naive parton model. In particular, in the one-photon approximation to deep inelastic 
lepton-hadron scattering, the relation ( |5.1|) is true to all orders in perturbation theory: 



N f 

F 2 (x,g 2 ) = £ e 2 \fT S) &Q 2 ) + fT s) &Q 2 ) ■ (5-31 



8=1 



Equivalently, one can say that in the DIS scheme the DIS coefficient functions are[| 

C? S{DIS \z;a s (Q 2 ) } Q 2 /ti = l) = C^ DIS \z;a s (Q 2 ),Q 2 /fi 2 F = l)=5(l-z) , 
Ci {DIS \z;a s (Q 2 ) } Q 2 /v 2 F = l) = 0. (5.32) 

In order to evaluate higher-order contributions in the small- a; regime, the DIS scheme 



offers some computational and phenomenological advantages [22]. The former amounts 
to the fact that in the DIS scheme one can explicitly resum the corrections as (as/N) k 
for the quark anomalous dimensions to all orders in as, as we shall illustrate below. 
As to the latter, we notice that the next-to-leading contributions as (as/N) k to the 
gluon anomalous dimensions are still unknown. Therefore, the knowledge of the quark 
anomalous dimensions in the DIS scheme may facilitate phenomenological investigations 
of the small- x behaviour of the structure function F2(x,Q 2 ). 

The comment above applies once the DIS scheme has been defined to all orders in 
perturbation theory. The point is that Eq. ( |5.31| ) (or, equivalently, ( p.32| )) fixes only the 
quark densities unambiguously. The relation between the singlet quark density in the DIS 
scheme and the MS-scheme parton densities is 

/gfV) = C 2 ^K(/i 2 ), 1) f s , N (fi 2 ) + C 9 2iN (a s (» 2 ), 1) iU(/i 2 ) (5.33) 

or, in terms of the matrix U a b in Eq. ( |5.26 ), 



Uss,N(a s ) = Cl N (a s A) , 2 N f U qg>N (a s ) = C{ N (a s , 1) . (5.34) 

The DIS-scheme gluon density, instead, still remains ambiguous and is given by an arbi- 
trary combination of gluon and singlet quark densities in the MS scheme 

V) = U gq , N (a s ^ 2 )) f s , N (v 2 ) + U mN {a s {^)) f g , N (» 2 ) , (5-35) 

with the only constraint ( |5.30|) , which reads 

U gg ,N=i(as) = 1 - C 9 2 N=l {a s , 1) , Ugq tN=1 (as) = 1 - Ci N (a s , 1) . (5.36) 



The convention introduced in Ref. [41] for defining fjj DIS ^ up to 0(as) amounts to 



extending Eq. ( 5.36|) to any value of iV in order as- A natural generalization of this 



convention is to require Eq. ( |5.36| ) to be valid for any N and to all orders in as- Doing 
that, the DIS-scheme gluon density is completely defined. 

Note, however, that for the purposes of our all-order calculation, it is not necessary 
to specify the actual form of the two matrix elements U gqi N and U g9t N in Eq. ( |5.35| ). We 



^Note that Eqs. ( 5.32 ) hold true only for a factorization scale fx 2 F = Q 2 
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just assume that they are chosen not to be extremely singular at high energies, i.e. they 
should not contain leading-order contributions of the type (as/N) k for iV — > 0. This is 
sufficient to ensure that most of the MS-scheme results obtained in the previous sections 
remain valid in the DIS scheme. In particular the gluon anomalous dimensions 7^ 75 ^ 

and the longitudinal coefficient functions C^ IS ^ are 



JPIS) 
lga,N 



C 



A (DIS) 
L, N 



-i ga ,N{a s ) + 0{a 2 s {a s /N) k 
Q 2 



g,q) 



0{a 2 s Q) k 



(A 



(5.37) 



where the resummed expressions for the MS-scheme anomalous dimensions 7 sa and co- 
efficient functions Cf are given in Eqs. (PH),(p^0|),(gllD, (|5^) . 



The quark anomalous dimensions, instead, do not coincide any longer (to this logarith- 
mic accuracy) with the corresponding anomalous dimensions in the MS scheme. Using 
Eqs. (|5~2TD and ( |533 ) we obtain 

1%*n\<xs) = J q9 , N (a s ) + <% s {ag, 1) 7 gg , N (a s ) + 0(a 2 s (a s /N) k ) . (5.38) 

On the other hand, the expression on the r.h.s. has been computed with logarithmic 
accuracy as(ois/N) k in Sect. 5.2. Inserting Eqs. (|5.17|) ,( |5.20| ) into ( |5.38|) , we find the 
following resummed expression for the quark anomalous dimension in the DIS scheme: 

2 + 3 7w - 3oi r'd - 7N ) m + 7w) RnM + o{o?s(as/Nf 



{DIS) i x Ot S „ 

l qg ,N (<x s ) = ^T d 
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3-2 lN r(2 + 2 77V )r(2-2 7 



N j 



The colour charge relation Q4.15 ) is still true in the DIS scheme: 



S(DIS) 
lqq,N 



(«s) 



C a 



7, 



(DIS) 
19, N 



(«s) 



as T 2 
2?r 3 



+ 0{a% (a s /N) k 



(5.39) 
(5.40) 



Some comments are in order. The DIS-scheme result (|5.39 ) for the quark anomalous 
dimensions has to be contrasted with the results discussed in Sect. 4 for the MS scheme. 
The algebraic complications of the MS scheme prevented us from obtaining resummed 
expressions in closed form for 7 g9 and Cf separately. We were able to explicitly resum 
only the combination in Eq. (|5.17|) , which turns out to be equivalent to the anomalous 
dimensions in the DIS scheme. This computational simplification has a more physical 
origin. The singlet sector of the deep inelastic lepton-hadron scattering is characterized 
by four physical observables, which can be studied in QCD perturbation theory: the 
structure functions F 2 S , Fj* and their first derivatives with respect to Q 2 . The MS 
scheme describes these observables in terms of eight different quantities: four coefficient 
functions Cf (i — 2, L, A = g, S) and the four matrix elements of the singlet anomalous 
dimensions. Obviously, only some linear combinations of them have to be regarded as 
physical observables. The DIS scheme, reducing to two the non-trivial coefficient func- 
tions (Cf = 1, Cf = 0), limits the number of arbitrary unphysical quantities necessary to 
describe the scattering process. The ensuing anomalous dimensions are more easily com- 
putable to all orders in as because they are more directly related to observable scaling 
violations. 



-33- 



Using the expansion ( 3.1 6| ) for the BFKL anomalous dimension 7 at and the expression 
(|3.17|) for R]y(as), the first perturbative terms of the quark anomalous dimension ( |5.39|) 
can be readily computed: 
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(5.41) 



The coefficients of the first two terms in the curly bracket agree with the one- and two- 
loop calculations in the DIS scheme [2^,28,30]. Moreover, all the coefficients in Eq. (|5.41|) 
are systematically larger than the corresponding MS-scheme coefficients in Eq. (|4.14j ). 
This behaviour is due to the additional contribution of the coefficient function Cf in 
Eq. (|5.38| ), and it is likely to persist in higher orders. 



Note also that the all-order expression ( |5.39| ) is analytic for < 7^ < 1/2. Thus, in- 
dependent of the value of as, the leading trajectory in iV-moment space is still given by 
the BFKL pomeron. As the BFKL anomalous dimension 7^ increases towards its satu- 
ration value at 77V = 1/2, the quark anomalous dimension quickly increases, approaching 
a singularity due to the pomeron normalization factor R^f(as) (see Eq. fl3.19|) ). This in- 
crease of 7^ 75 ^ leads to strong scaling violations, although the singularity at 7^ = 1/2 is 
cancelled in physical observables by analogous contributions to the resummed coefficient 
functions in Eq. ( [5. 37] ). 



6. Summary 

In the present paper we have shown how the high-energy factorization theorem H 
can be extended beyond the leading logarithmic accuracy in a manner which is consistent 
with the all-order factorization of collinear singularities. Much effort has been devoted to 
investigating the issue of the dependence on the factorization scheme of parton densities 
and coefficient functions. This analysis has led to the (off-shell) fc^-factorization in di- 
mensional regularization represented schematically by Eq. ( |2.29| ) (see also Eqs. ( [4.2]) and 
~])). Eq. ( E39D has then been used to study the high-energy (or small-x) behaviour of 



deep inelastic scattering processes. 

A first general consequence of Eq. Q2.29|) is that flavour non-singlet observables are 
regular at small x order by order in perturbation theory. 



As regards the singlet sector, fcj_-factorization allows one to sum classes of logarithmic 
corrections to all orders in as- To this end, one has to evaluate 2GI kernels (see Eq. (|2.29|) ) 
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in fixed-order perturbation theory and use the master equations (|3~2"|),( |3~!?(]| ) for the gluon 
Green functions. 



Eqs. ( |3.2| ) and ( |3.20| ) are the generalization of the BFKL equation [T2|] to the case 
of n = 4 + 2 e space-time dimensions. Their solution is discussed in Sect. 3. In par- 
ticular, the calculation of the gluon anomalous dimensions 7 95i jv(a:s), 7 ot ,at(q;s) (see 
Eqs. ( ft.l3D ,( |O4|) ,( |3.30D ) has been carried out to the leading logarithmic accuracy (as/N) k 
in the context of dimensional regularization and the BFKL results on the pomeron tra- 
jectory have been re-derived. Besides this, we have been able to compute (to the same 
accuracy) the normalization factor R^^as) (see Eq. ( [3.17f) ) of the perturbative QCD 
pomeron in the MS factorization scheme. 

The quark sector enters the QCD evolution equations to next-to-leading logarithmic 
order as {as lnx) fc . Using high-energy factorization, in Sect. 4 we have evaluated the cor- 
responding quark Green functions. We have shown that the result in Eq. (|4.1U[ ) originates 
from an integral equation whose kernel is related to a generalized (off-shell) Altarelli- 
Parisi splitting function (see Eq. ( f4.9|) ). We have also discussed how Eq. ( |4.10|) can be 
used for evaluating the small- N limit of the quark anomalous dimensions 7 g9 , jv(a;,s*) and 
1qq,N( a s)- The result of our explicit calculation up to six- loop order is given in Eq. (|4.14j ). 

As an example of application of high-energy factorization to a specific hard process, 
in Sect. 5 we have considered deep inelastic lepton-hadron scattering (for the case of 
heavy-flavour production see Refs. and Sect. 2.4). Resummed expressions to next-to- 
leading accuracy as(®.s/N) k for the DIS coefficient functions C 2 and Cx are presented 
in Eqs. (|5TTTD , (gT7|) , (|5^I|) , (gjg) , (^23|) . These results are given in the MS factorization 



scheme. In Sect. 5.3, we have also introduced an all-order generalization of the DIS 
factorization scheme first proposed in Ref. f4T| . Within this scheme, where most of the 
DIS coefficient functions are trivial (see Eqs. (|5.32|) and ( |5.37|) ), we have obtained the 



next-to- leading resummed expressions ( |5.39| ),( |5.40| ) for the quark anomalous dimensions. 

Quantifying precisely the phenomenological consequences of the results presented here 
is a matter of detailed numerical investigations. However, the size of the next-to-leading- 
order coefficients in the perturbative expansions (|4.14j ),( |5.24P , (|5.25|) and (|5.41| ) suggests 



that these contributions may have phenomenological relevance in accurate analyses of scal- 
ing violations, already at the values of a; (x ~ 10~ 3 -r- 10 -4 ) accessible at present hadron 
colliders. In particular, since the first leading-order coefficients of the gluon anomalous 
dimensions are vanishing (see Eqs. (|3.16|) ), the next-to- leading corrections in the quark 



sector computed in this paper may be quite important for the study of the proton struc- 
ture functions being measured at HERA. A fully consistent analysis to next-to-leading 
logarithmic order obviously requires also the computation of the still unknown (to this 
accuracy) gluon anomalous dimensions. We hope to report progress on this subject in the 
near future. 
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Appendix A 



The rescaled parton densities f a (x,fi 2 ) = x f a (x, /z 2 ) introduced in Eq. ( |2.6| ) fulfil the 
evolution equations 



d fafafJ? 



V/ dz P ab (a s (fi 2 ),z) f b (x/z,fi 2 



(A.l) 



where P a b(as,z) are the generalized Altarelli-Parisi splitting functions in Eq. ( |2.3| ). Ac- 
cording to our normalization, the splitting functions have the following one-loop expres- 
sions 



2C A 
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= T R 
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1 + + z(l-z) 
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(A.2) 



in terms of the SU(N C ) colour factors (N c = 3 is the number of colours) 



C A =N n 



N 2 c -l 

2 N r . 



Tr (t a t b )=6 ab T R = -6 ab 



(A3) 



The leading-order splitting functions P^ b '(z) (which have been known for a long time 
||26|| ) are factorization theorem invariants, i.e. they do not depend on the explicit procedure 
to factorize collinear singularities. The physical reason for this is that they are directly 
related to observable scaling violations in deep inelastic scattering processes. On the 
contrary, splitting functions and anomalous dimensions beyond one-loop order do depend 
on the regularization and factorization schemes of collinear singularities. Nonetheless, due 
to charge conjugation invariance and SU(Nf) flavour symmetry of QCD, they satisfy the 
following scheme-independent properties 

lq%9 = ^qt9 = Iqg > Igm = 1gq~i = Igq 

Tfcfc = iMj = IgY^j + 7« , lm = lm = Iq^H + Iqq • ( A - 4 ) 



The symmetry properties ( |A.4[ ) imply that the anomalous dimensions matrix has 
only seven independent components. Correspondingly, three flavour non-singlet (f^ v \ 
fq~, fqf) an d two flavour singlet (fs, f g ) parton densities can be introduced so that the 
evolution equations ( |A.1| )) (or ( |2.9| )) are completely diagonalized (in the partonic space) 
for the non-singlet sector. One explicitly finds (we drop the overall dependence on N, // 
and as) 

-j\ — r = T ' i] % = ^ ' i\ 2 = ^ ' ( A - 5 ) 

a in /i z a In [x A a in 
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where 



Nf 1 Nf 

?(v) = V7f _ f _ ^ f(±) _ » ■• 1 



/ — /<2j) ) /oj — /?i ^ /?i AT y^Xfcin ^ 5 (^--6) 

i=l ^ 3=1 

and the non-singlet anomalous dimensions are given by 

1 {V) = li S - 7 g f + N f (^ g - 7 * ) , 7 W = 7? f ± 7? f . (A.7) 
The evolution equations are instead still coupled in the singlet sector: 
dfs 



dln/[i 2 " K S + ^ S + N f^ + &s + 2N flqg f g 

= ljs + lj 9 , (A. 



where the quark singlet density is defined by fs = YnJiifqi + 

From Eq. ( |A.2| ) we see that all the three non-singlet anomalous dimensions are de- 
generate in one-loop order. The anomalous dimensions in two-loop order were computed 
in Refs. In this order the degeneracy mentioned above is partially removed 

because 7 ^' s ^ 0. However we still have 7^'(ag) = 7^(0:3) + 0(a 3 s ) since 7 ^ g and 7 ^ 
coincide in 0{a 2 s ). The equality between 7 ^ g and 7 ^ is expected to be violated starting 
from 0(a\). 

The high-energy power counting in Sect. 2.3 implies that the non-singlet anomalous 
dimensions in Eq. ( |A.7| ) (and hence 7^jvj ^qq S Ni IqqN ~ 7gg,Ar) are regular for A" — > 
order by order in as- In other words, the corresponding n-loop splitting functions p( n ~ x ) 
in Eq. (|2.3| ) are less singular than 1/x for x — * 0. All the high-energy contributions 
Og/N k {n > k > 1) are thus associated with the gluon anomalous dimensions •y gg ,N,lgq,N 
and with the quark anomalous dimensions 1q g ,N,1 ggj N — 1qq,N- 
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Appendix B 



In order to obtain the solution of the master equation ( |3.2| ) for the /c^-distribution 
it is convenient to introduce the dimensionless distribution T. as follows 



F$\k;a s ,»,e) = 5^\k) + fA±Z±A f N (a s (k 2 /^ 2 y,e) . 



ra + e) 



7T k 

Note that f does not depend on e, as and k independently. In fact, rewriting Eq. fl3T2] 
in terms of JF, we see that it only depends on e and the combination as(k 2 / fx 2 ) £ : 
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(B.2) 



Then, we perform the shift q — > q' = q — k of the integration variable in Eq. ( |B.2| ). 
and introduce the ratio r = |q'|/|k| and the angle 9 between q' and k. The integral 
equation ( p.2| ) now reads: 
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The angular integration in Eq. ( |B.3| ) can be performed in terms of associated Legendre 
functions. However we did not find the ensuing representation of any convenience. 



The perturbative solution of Eq. 

F N (a s (k 2 /v 2 y,e) = Y: 



k=l 



has the following expression: 

'l.2\ 



«5 

n e rn + £ 



Inserting Eq. (|B.4[ ) into Eq. ( |B.3| ), we obtain the recurrence relation 
bative coefficients Cfc(e). The recurrence factor h{s) is given by 

I k (e) = 7(7 = ke;e) , 

where J( 7 ; e) is the following integral 
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(B.6) 
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Using Eqs. (B.5) and 



we obtain Eq. 



represents the action of the kernel of the master 

2\ 



The integral 1(7; e) in Eq. (|B 
equation ( |B.2| ) (or, equivalently, (|3.2|)) on a test function .F(k) behaving like (k z ) 7 . 
Therefore the fact that 1(7; e) is finite for e — » and 7 fixed is a consequence of the 
collinear regularity of the kernel of the BFKL equation in four dimensions. Note also 
that this regularity is achieved through the cancellation of collinear singularities which 
are present separately in the real and virtual contributions, i.e. in the first and second 
terms in the square bracket of Eq. ( |B.2| ) or (p.6| ) . 



As discussed in Sect. 3, the lack of scale invariance in the master equation ( |B.2| ) does 
not allow us to find a general solution for T for any value of e. This means that we are 
not able to resum the formal power series expansion in Eq. ( |B.4| ). Nevertheless, we can 
obtain an explicit all-order solution in the relevant limit of small e values. This limit is 
sufficient to compute the transition function T gg (i.e. the gluon anomalous dimensions) 
and the associated normalization factor in Eqs. ( p.8|) ,(|3~9|). 



To do this, let us perform the derivative of Eq. 
divide both sides by T . We thus obtain 
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Unlike Eq. ( |B.3| ), Eq. ( |B.7| ) is homogenous with respect to T . Therefore we can easily 
factorize the singular transition function Yg g ^{as-, e). 

More precisely, let us notice that the distribution T introduced in Eq. ( |B.1| ) is related 
to the integrated gluon Green function (|3.6| ) as follows 
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Using the factorization formula 
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we thus have 
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Since the anomalous dimensions 7j\f(«s 5 E ) and the renormalized Green function 
G gg , n(ois, e) are regular for e — > 0, it follows from Eq. ( |B.10|) that the relevant func- 
tions involved in Eq. ( |B.7|) and Eq. (|B.7l ) itself can be expanded in e around e = 0. A 
straightforward calculation gives: 
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where Rn(qs) = G ggj iy(as,£ = 0). Correspondingly, Eq. ( |B.7| ) reads as follows (we drop 
the overall explicit dependence on as) 
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where we have introduced the (second-order) differential operator Q 
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acting on the integral I^^e) in Eq. ( |B.6| ). Equating the O(e ) and 0{e) terms in 
Eq. ( |B.12|) we thus find an implicit equation for 7^(0:5) and a differential equation for 
R N (a s ). 



To order e° , we obtain 
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i.e. the result (|3.14 ) for the BFKL anomalous dimension, x{l) being the characteristic 
function in Eq. ( |3.15 ). To order e , we have 
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where x'(7) > x"{l) are the first and second derivatives of x(l) with respect to 7. In 
order to solve Eq. ( [B.15|) , it is convenient to consider 7^ as the independent variable 
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and use Eq. (|B.14[) : 
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Inserting Eqs. (p.l6| ),( [B.17|) into Eq. (p.!5| ) we get the differential equation 
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whose solution is given by Eq. ( |3.17|) . 



Appendix C 



The evaluation of the quark anomalous dimensions starting from the fcj_-factorization 
formula ( |4.4|) requires the explicit computation of the off-shell kernel K qg . Inserting 
Eq. ( |4.2| ) into Eq. ( |4.3| ) and comparing the latter with Eq. ( |4.4j ) we find 
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where K qg (z) is obtained from K^\q, k) (see Fig. 7b), after integration over q, as follows 
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In Eq. ( |C.2|) the following Sudakov parametrization for the momenta k and q is un- 
derstood: 

2 2 

k» = yp» + k t \ , q>* = xp l * + q>t + l , z=- . (C.3) 

2 x p ■ p y 

It is also convenient to introduce the boost-invariant (along the fc-direction) transverse 
momentum q 

q = q-^k . (C.4) 
Performing the Dirac and colour algebra we get 
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and, inserting Eq. ( U.5 ) into Eq. ( |C2| ), the azimuthal average over q and the integration 
over q 2 can easily be performed, thus leading to 
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Eq. 



is precisely the result in Eq. 



expressed in terms of the off-shell splitting 



function 



The q 2 -integration in Eq. ( |C.6| ) is not elementary and provides a representation of the 
off-shell kernel K qg in terms of hypergeometric functions. However, in order to obtain 
the power series expansion ( |4.10| ) for the Green function G q ° g \ it is more convenient to 
carry out first the k-integration in Eq. (|4.4|) . Inserting the expansion ( |3.3|) into Eq. (|4.4j ), 

we get 
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where the function /i gfl ,jv(7;e) is defined by the following /c^-transform of the off-shell 
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The evaluation of /i gg is straightforward. We insert Eq. ( p.6| ) into Eq. ( |C.8| ) and perform 
first the trivial integration over k 2 with q 2 /(,z(l — z)~k 2 ) = r fixed. Then, the integrations 
over t and z decouple and can be carried out in terms of Euler gamma functions. Note 
also that the N — > limit of h qg ^ is regular (i.e., K qg (z) ~ z, modulo lnz-terms, for 
z — > and any value of k 2 /Q 2 ) . Therefore, to the logarithmic accuracy we are interested 
in, we can limit ourselves to computing h q9j N=o- The final result is: 
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Using Eq. ( |C9|) , we recover Eqs. ( 4.10 ) and ( [4.11 ) by the identification 
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As discussed in Sect. 4, we have not been able to use the power series expansion 
(|4.10 ) for explicitly resumming all the next-to-leading logarithmic corrections in the quark 
anomalous dimensions r y qg ,N- in general, to this accuracy we can write 
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and the values of the coefficients a k for k < 5 are given in Eq. Q4.14| ). The calculation of 
the higher-order coefficients is much more cumbersome. Here, we present only the result 
for the rational part of a k . In other words, let us split a k as follows 



ak = r k + b k 



(C.12) 



where b k is an irrational number given in terms of powers of Riemann zeta functions 
{n > 3) and r k is the residual rational contribution to a k . We find the following expression 
for r k in the MS scheme: 
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One can easily check that Eq. ( |C.13| ) reproduces the values of r k in Eq. ( f4.14| ), i.e. r\ = 
5/3, r 2 = 14/9, r 3 = 82/81, r 4 = 122/243, r 5 = 146/729. The derivation of the result 
(|C.13| ) is left as an exercise for the reader. 
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Figure Captions 



Figure 1: Factorization of the physical cross section F in terms of partonic cross sections 
F 4 (0) , . . .) and parton distributions (/(°>, fi°\ ...). In the case of on- 
shell partons (pf = 0), the first term on the r.h.s. and those in parenthesis 
represent respectively the leading- and higher-twist contributions. 

Figure 2: Expansion of the partonic cross section F^ in two-particle irreducible (2PI) 
kernels. 

Figure 3: Action of the collinear projector Vc on the spin indices of two kernels A and 
B in the cases of an intermediate (a) quark and (b) gluon state. 

Figure 4: Expansion in two-gluon irreducible (2GI) kernels at high energy for (a) the 
partonic cross section F^> and (b) the (singlet) quark Green functions Q qa , 

Qqa- 

Figure 5: Quark-antiquark contribution to the lowest-order absorptive part A^ u of the 
scattering amplitude 7 g — > 7 g. 

Figure 6: The BFKL characteristic function ^(7) for < 7 < 1. 7^ is the BFKL 
anomalous dimension. 

Figure 7: (a) High-energy factorization of the gluon — > quark Green function Q qg and 
(b) the corresponding off-shell kernel K qg . 
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